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Summary 

Even  in  the  absence  of  point-and-tracking  turrets,  airborne  lasers  still  need  to  propagate 
through  turbulent  boundary  layers  present  on  a  skin  of  an  aircraft.  Previous  AFOSR- 
funded  work  (Grant  FA9550-09- 1-0449)  performed  systematic  studies  of  the  optical 
character  of  attached  turbulent  boundary  layers.  These  studies  also  included  non-adiabatic, 
primarily  heated  walls. 

This  prior  experimental  work  exploring  the  aero-optics  of  boundary  layers  has  shown 
strong  evidence  for  the  dominance  of  large-scale  turbulent  structures  as  a  main  source  for 
boundary  layer  induced  aberrations.  So,  one  way  to  mitigate  aero-optical  effects  is  to 
reduce  the  size  of  the  large-scale  structures.  One  of  the  devices  successfully  used  to 
passively  manipulate  the  large-scale  structure  is  a  Large-Eddy  Break-Up  (LEBU)  device. 
LEBU  devices  were  shown  to  be  effective  at  reducing  streamwise  and  spanwise  integral 
length  scales,  boundary  layer  intermittency,  and  the  strength  of  bulges  and  sweeps  without 
introducing  a  strong  shear  layer  or  other  optically  un-desirable  flow  features,  making  them 
an  attractive  candidate  for  use  as  an  aero-optic  mitigation  application  for  TBLs.  While  the 
potential  for  using  LEBUs  for  aero-optic  mitigation  was  recognized  over  three  decades  ago 
these  suggestions  have  not  been  seriously  investigated. 

Another  approach  to  mitigate  aero-optical  effects  is  to  reduce  the  density  fluctuations. 
In  the  previous  work  it  was  also  shown  that  by  heating/cooling  the  wall  the  density 
structure  can  be  manipulated  without  any  significant  effects  on  the  velocity  structure. 

In  the  present  work,  the  effect  of  two  passive  boundary  layer  flow  control  techniques, 
LEBU  devices  and  wall  heating/cooling,  on  turbulent  boundary-layer-induced  aero-optical 
aberrations  is  experimentally  investigated.  Extensive  parametric  studies  are  performed 
investigating  the  effect  of  LEBU  geometry  on  levels  of  optical  aberrations  in  the  turbulent 
boundary  layer.  The  results  of  these  experiments  are  analyzed  to  determine  the  physical 
mechanisms  responsible  for  the  experimentally  observed  changes.  It  was  demonstrated  that 
long  LEBUs,  both  single  and  multi-element  ones,  are  effective  passive  means  to  reduce 
aero-optical  distortions  by  as  much  as  45%  immediately  downstream  of  the  device,  with 
the  effect  lasting  for  several  boundary  layer  thicknesses. 

The  effect  of  moderate  levels  of  boundary  layer  wall  cooling,  both  for  full  and  partial 
wall  cooling,  on  aero-optic  aberrations  is  also  experimentally  investigated.  The  results  are 
favorably  compared  to  a  statistical  model  derived  using  the  temperature-velocity  relation 
from  the  Extended  Strong  Reynolds  Analogy.  For  the  optimal  cooling  the  aero-optical 
effects  were  demonstrated  to  be  reduced  by  60%.  It  was  shown  that  the  cooling  might  be 
also  effective  in  reducing  aero-optical  distortions  of  the  boundary  layers  at  supersonic 
speeds.  Comparing  the  model  prediction  to  the  experimental  results,  it  was  speculated  that 
while  the  pressure  effects  can  be  neglected  in  adiabatic  boundary  layers,  they  might  play 
an  important  role  in  moderately-cooled  boundary  layers. 

As  moderate  heating  simply  introduces  passive  temperature  markers  in  the  boundary 
layer  and  effectively  amplify  aero-optical  levels,  by  introducing  moderate-temperature 
mismatch  between  the  wall  and  the  boundary  layer,  it  is  possible  to  thermally  tag  and 
measure  three-dimensional  large-scale  structures  in  incompressible  boundary  layers.  This 
important  result  allows  using  non-intrusive  wavefront  sensors  as  alternative  complimentary 
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diagnostic  tool  and  in  the  present  work  it  was  used  to  study  aero-optical  properties  of 
turbulent  boundary  layers  with  low  Reynolds  numbers. 

Finally,  simultaneous  velocity/wavefront  measurements  were  performed  in  both 
incompressible/heated  and  compressible  boundary  layers  in  order  to  demonstrate  that 
optical  sensors,  combined  with  traditional  measurement  techniques,  can  be  useful  in 
studying  fundamental  aspects  of  physics  and  dynamics  of  large-scale  structures  in 
boundary  layers.  Analysis  of  instantaneous  velocity  and  density  fields  leads  to  a  better 
understanding  of  the  connection  between  aero-optical  distortion  in  turbulent  flows  and 
turbulent  velocity  structure.  By  comparing  the  velocity  and  optical  data  and  using  an 
instantaneous  version  of  Strong  Reynolds  Analogy,  plausible  relations  between  different 
large-scale  structures  in  the  boundary  layer  and  the  corresponding  wavefronts,  as  well  as 
the  effect  of  the  local  pressure  field  inside  the  vortical  structures,  are  presented  and 
discussed. 
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1.  Aero-Optics  of  Boundary  Layers 


Within  the  aerospace  community,  there  is  substantial  interest  in  developing  aircraft-based 
directed  energy  systems  for  a  variety  of  applications,  including  air-to-air  and  air-to-ground,  line- 
of-sight  optical  telecommunications  systems.  This  type  of  system  would  use  laser  beams  to  relay 
data  directly  between  optical  transmitters  and  receivers  on  aircraft,  ground  stations,  and  satellites, 
and  would  allow  for  much  greater  communication  bandwidth  than  is  currently  available  on  state- 
of-the-art  aviation  communication  technologies. 

One  of  the  many  challenges  that  must  be  overcome  before  line-of-sight  optical 
communications  systems  can  become  a  robust  and  reliable  method  of  communication  for  airborne 
vehicles  is  the  effect  that  the  Earth’s  atmosphere  has  on  the  propagation  of  light.  This  fluid-optic 
interaction  occurs  when  wavefronts  propagating  through  the  Earth’s  atmosphere  encounter 
regions  of  air  that  have  non-uniform  and/or  unsteady  distributions  of  density.  Since  density  and 
index-of-reffaction  are  proportional  to  one  another,  these  regions  of  fluctuating  density  cause 
initially  planar  optical  wavefronts  passing  through  them  to  be  distorted  (Gladstone  &  Dale,  1863; 
Liepmann,  1952;  Tatarski,  1961,  1971).  These  distortions  can  cause  wavefronts  to  destructively 
interfere  with  themselves,  leading  to  significant  reductions  in  the  average  on-target  beam  intensity 
that  reduce  the  overall  performance  of  directed  energy  systems  (Gilbert,  2013;  Bom  &  Wolf, 
1999;  Jumper  &  Fitzgerald,  2001).  Unsteadiness  in  the  density  field  has  also  been  linked  to 
instantaneous  reductions  in  beam  intensity,  which  can  cause  disastrous  signal  drop-outs  for  free- 
space  communications  applications  (Gordeyev,  Cress,  and  Jumper,  2013).  Wavefront  distortions 
that  are  the  result  of  density  fluctuations  found  in  regions  of  turbulent  flow  surrounding  aircraft 
are  identified  as  the  aero-optic  problem  (Gilbert  &  Otten,  1980;  Jumper  &  Fitzgerald,  2001; 
Wang,  Mani  &  Gordeyev,  2012). 

The  levels  of  optical  wavefront  distortions  can  be  quantified  by  the  Optical-Path- 
Difference,  OPD(x,z,t),  (Wang  at  el,  2012) 

b 

OPL(x,  y,t)  =  KGD\  p'  (x,  y,  z,  t)dy 

J  n  n 


OPD(x,  y,  t )  =  OPL(x,  y,  t )  -  ( OPL{x ,  y,  t)) 


where  KGD  is  the  Gladstone-Dale  constant,  the  integration  is  performed  along  the  beam 
propagation  axis,  y,  and  the  angled  brackets  denote  spatial  averaging.  The  extent  to  which 
wavefronts  are  distorted  from  their  initially  planar  state  is  typically  characterized  by  taking  the 
time-averaged  root-mean-square  (RMS)  of  the  OPD,  or  OPDrms. 

The  performance  of  free  space  optical  systems  can  be  meaningfully  quantified  by  defining 
the  Strehl  ratio,  SR  =  I/I0,  where  I  is  the  on-target  beam  intensity  in  the  far  field,  and  Io  is  the 
diffraction-limited  performance  in  the  far  field.  In  the  presence  of  optical  distortions  SR<1,  and, 
for  small  levels  of  OPDrms  for  a  given  laser  wavelength,  X,  can  be  estimated  as, 


SR  =  exp 


2  nOPD^  \ 

2 

K  X  j 

so-called  Large- Aperture  Approximation  (Ross,  2009;  Porter,  et  al.  2013) 

Sutton  (1969)  introduced  the  most  widely  cited  theoretical  formulation  for  calculating  the 
effect  of  turbulent  boundary  layers  on  optical  wavefronts  from  turbulence  statistics.  The  equation 
derived  by  Sutton,  referred  to  as  the  aero-optic  ‘linking  equation,’  was  based  heavily  on 
Tatarski ’s  (1961)  work  describing  the  distortions  of  electromagnetic  waves  propagated  through 
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the  Earth’s  atmosphere.  It  related  the  time-averaged  levels  of  OPDrms  to  turbulence  density 
statistics.  A  simplified  version  of  the  linking  equation  is 

OPDL=2K2eDjpl(y)Al,(y)dy,  (1.2) 

where  p2ms  is  the  variance  of  the  density  fluctuations,  and  A  p  is  the  local  density  correlation 

length.  The  linking  equation  (1.2)  has  been  validated  both  experimentally  (Hugo  &  Jumper,  2000; 
Gordeyev  et  al,  2014)  and  numerically  (Wang  &  Wang,  2013;  Wang  &  Wang,  2013)  and 
establishes  a  relation  between  the  level  of  aero-optical  distortions  and  properties  of  the  density 
structure,  /7nils  and  Ap.  So,  aero-optical  data  can  provide  important  information  about  the  density 
fields  in  turbulent  flows. 

In  order  to  achieve  a  large  field  of  regard,  airborne  directed  energy  systems  have 
gravitated  towards  the  use  of  hemisphere  or  hemisphere-on-cylinder  turrets  in  order  to  send  and 
receive  laser  beams.  While  these  devices  are  mechanically  simple,  the  complex  three-dimensional 
flow  that  develops  around  turrets  includes  elements  of  turbulent  boundary  layers,  flow  separation, 
shear  layers,  wakes,  and  unsteady  shocks  (Gordeyev  &  Jumper,  2010). 

For  applications  in  which  a  large  field  of  regard  is  not  an  essential  requirement,  the  flow- 
field  over  the  optical  aperture  can  be  vastly  simplified  by  eliminating  the  turret  completely. 
Instead,  the  beam  could  be  propagated  through  an  aperture  or  phased  array  that  is  flush-mounted 
to  the  skin  of  the  aircraft  (Serati  &  Stockley,  2003;  Whiteley  &  Gordeyev,  2013).  In  this 
configuration,  the  only  aero-optically  active  flow  field  that  must  be  traversed  by  the  beam  is  the 
compressible,  turbulent  boundary  layer  (TBL).  At  present,  aero-optic  aberrations  caused  by  the 
boundary  layer  have  been  studied  analytically,  experimentally  and  computationally  for  over  six 
decades  (Gordeyev  et  al.  2014  and  references  therein),  and  much  has  been  learned  about  the 
physics  of  TBL  fluid-optic  interactions.  Important  findings  from  research  on  the  aero-optical 
characteristics  of  TBLs  can  be  summarized  via  a  few  main  conclusions.  First,  aero-optic 
aberrations  have  been  strongly  correlated  with  the  large-scale  turbulent  motions  convecting  in  the 
outer  part  of  the  TBL  (Gordeyev  et  al.  2014,  and  the  references  therein).  Secondly,  analysis  of 
time-resolved  TBL  wavefront  measurements  has  revealed  that  even  in  cases  where  on-average 
target  intensity  is  found  to  be  high,  large  intermittent  drop-outs  in  far-field  signal  intensity  that 
lasted  on  the  order  of  one  millisecond  were  found  (Gordeyev,  Cress  &  Jumper,  2013),  which 
would  be  crippling  to  the  performance  of  any  potential  airborne  laser-based  high-bandwidth  free- 
space  communications  systems. 

These  results  clearly  demonstrated  that  although  the  TBL  is  vastly  simpler  than  the  flow 
environment  around  a  turret  body,  boundary  layers  are  not  guaranteed  to  be  entirely  benign  with 
respect  to  aero-optics.  As  a  result,  there  is  substantial  interest  in  developing  and  testing  flow 
control  schemes  that  would  modify  the  boundary  layer  flow  in  such  a  way  that  would  permit 
airborne  directed  energy  and  free-space  communications  systems  to  achieve  high  on-target 
intensity  in  both  a  time-averaged  and  instantaneous  sense. 

The  linking  equation  (1.2)  implies  that  in  order  to  reduce  overall  aero-optical  distortions 
of  a  boundary  layer,  we  should  investigate  flow-control  approaches  that: 

1 .  Reduce  density  fluctuations  inside  the  boundary  layer,  primarily  in  its  outer  region, 

2.  Reduce  the  size  of  aero-optical  structure, 

3.  Reduce  the  integration  region. 
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Jumper  (1980)  first  suggested  the  use  of  parallel  plate  manipulators  for  aero-optic  flow 
control.  These  devices,  which  are  also  commonly  referred  to  as  Large-Eddy  Break-Up  (LEBU) 
devices,  were  originally  proposed  as  a  method  of  flow  control  for  achieving  viscous  drag 
reduction  (Corke,  et  al.  1979).  LEBUs  were  found  to  achieve  this  end  by  altering  the  turbulence 
structure  in  the  outer  part  of  the  TBL,  and  they  significantly  reduced  integral  length  scales  for 
long  distances  downstream  of  the  devices  (Corke,  1981;  Savill  &  Mumford,  1988;  Anders,  1990, 
and  others).  These  effects,  especially  the  suppression  of  large-scale  turbulent  structures  in  the 
outer  part  of  the  boundary  layer,  are  very  much  aligned  with  the  mechanisms  identified  as  the 
dominant  source  of  TBL  aero-optic  aberrations.  Prior  to  this  work  however,  no  experimental  or 
computational  evaluation  of  LEBUs  for  aero-optic  mitigation  has  been  performed.  Chapter  2  of 
this  report  will  cover  main  results  of  this  study. 

In  fact,  very  few  studies  of  flow  control  methods  for  mitigating  TBL  aero-optics  have 
been  performed  to  date  (Cress,  2010;  White  &  Visbal,  2011).  One  illuminating  investigation  was 
performed  by  Cress  (2010)  under  AFORS-funded  Grant  FA9550-09- 1-0449,  in  which  the  effect 
of  wall-temperature  on  TBL  wavefront  aberrations  was  experimentally  investigated.  Experiments 
were  particularly  focused  on  wall-heating,  but  both  preliminary  experimental  measurements  and  a 
simplified  theoretical  model  for  wall  temperature  effects  indicated  that  wavefront  aberrations 
could  be  reduced  by  as  much  as  80  %  through  the  use  of  wall  cooling  (Cress,  2010).  Additional 
data  are  still  required,  however,  to  validate  and  improve  Cress’  model  over  a  wide  range  of 
temperatures  and  Mach  numbers,  but  these  findings  already  suggested  that  wall  cooling  is  also  a 
promising  technique  for  aero-optic  mitigation.  Extensive  work,  including  a  development  of  the 
theoretical  model  extended  to  supersonic  regime  was  done  during  this  grant  and  overview  of 
major  results  is  presented  in  Chapter  3. 

Also,  there  is  a  pressing  need  to  develop  techniques  to  ‘close  the  gap’  in  Reynolds  number 
between  experiments  and  computations  of  the  aero-optic  effects  of  TBLs.  Due  to  both 
computational  constrains  and  wavefront  sensor  limitations,  experimental  measurements  of  TBL 
aero-optic  effects  are  performed  at  Reynolds  numbers  that  are  about  an  order  of  magnitude  higher 
than  results  from  high-fidelity  simulations,  making  one-to-one  validation  of  the  computational 
methods  used  difficult  to  date.  Wall  heating  experiments,  presented  in  details  in  Chapter  4, 
showed  that  it  is  possible  to  passively  amplify  TBL  wavefront  aberrations,  thus  raising  the 
possibility  of  using  wall  heating  as  a  means  of  increasing  the  signal-to-noise  ratio  in  low-speed 
experiments  where  aero-optic  aberrations  are  otherwise  very  weak. 

Finally,  non-intrusive  feature  of  optical  measurements,  as  well  as  sensitivity  to  density 
field  only,  makes  them  very  attractive  complimentary  diagnostic  tools.  When  combined  with 
traditional  measurement  techniques  like  hot-wires,  PIV,  pressure  sensors  etc,  they  provide 
valuable  additional  information  about  physical  mechanisms  inside  the  boundary  layers.  To 
demonstrate  that,  simultaneous  velocity/optical  measurements  were  performed  in  both  low 
subsonic  and  transonic  boundary  layers.  By  analyzing  corresponding  velocity  and  density  fields 
and  assuming  different  relations  between  these  quantities,  a  better  understanding  of  the  physics 
and  dynamics  of  the  large-scale  vortical  structures  and  their  role  in  aero-optical  distortions  can  be 
obtained.  In  particular,  preliminary  results  indicate  that  the  pressure  field  might  not  be  negligible 
inside  these  structures.  A  summary  of  these  results  are  provided  and  discussed  in  Chapter  5. 
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2.  Aero-Optical  Mitigation  Using  Large-Eddy  Break  Up  Devices 


2.1  Motivation 


A  great  deal  of  effort  in  the  field  of  fluid  mechanics  has  been  dedicated  to  identifying  and 
evaluating  different  methods  of  passive  flow  control  for  turbulent  boundary  layers,  including 
riblets  (Walsh  1990),  vortex  generators  (Zaman,  Hirt,  &  Bencic,  2012),  honeycombs  and  screens 
(Yajnik  &  Ancharya,  1977),  and  large-eddy  break-up  (LEBU)  devices  (Corke,  Guezennec,  & 
Nagib,  1979).  Studies  of  LEBU  devices  for  reducing  viscous  drag  in  turbulent  boundary  layers 
has  confirmed  that  C/  could  be  reduced  for  distances  -O  (100b)  downstream  (Corke,  1981; 
Anders,  1990)  by  substantially  altering  the  turbulence  structure  in  the  outer  part  of  the  boundary 
layer  (Anders,  1990).  The  effectiveness  of  the  LEBU  device  for  modifying  large-scale  turbulent 
structures  in  the  outer  part  of  the  TBL  has  caused  it  to  be  identified  as  a  potential  method  of 
passive  flow  control  that  could  be  used  to  achieve  aero-optic  mitigation  in  the  TBL  (Jumper, 
1980;  Anders,  1985). 

In  general,  LEBU  devices  consist  of  one  or  more  thin  flat  plates  or  airfoils  placed  parallel 
to  the  wall  within  the  turbulent  boundary  layer,  as  shown  in  Figure  2.1.  These  devices  were 
designed  to  act  as  screens  with  sparse  elements,  significantly  reducing  blockage  effects  while 
retaining  the  ability  to  directly  manipulate  turbulent  structures  in  the  outer  part  of  the  boundary 
layer  (Corke,  1981).  The  significant  reduction  in  blockage  compared  to  screen-type  devices  was 
expected  to  result  in  a  net  viscous  drag  reduction  compared  to  the  un-manipulated  TBL  (Corke, 
Guezennec,  &  Nagib,  1979).  Some  experiments  have  investigated  LEBU  devices  with  as  many  as 
four  flat  plate  elements  (Corke,  1981),  however  this  work  is  limited  to  investigating  LEBU 
devices  with  at  most  two  flat  plates. 


Cf  reductions  persist 
for  0(1 00r)')  r— ^ 


Tapered  Leading  &  Trailing  Edges 


LEBU  Device 


.  LEBU- 
v  Modified  TBL 


Device  wake 


Incoming  Turbulent  Boundary  Layer 

m 


Flow 


Flow  visualization  from  Corke,  et  al.  (1979) 


a)  Single  LEBU  b)  Multi  LEBU  c)  Tandem  LEBU 


Figure  2.1.  Top:  Schematic  of  the  effect  of  a  single  LEBU  device  on  the  TBL.  Bottom: 
schematics  of  the  a)  Single  LEBU,  b)  Multi-LEBU,  and  c)  Tandem-LEBU  configurations,  [flow 
visualization  image  from  Corke,  et  al.  (1979)]. 
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Single  LEBU  devices,  shown  in  Figure  2.1a,  consist  of  a  single  thin  plates  or  airfoils  of 
thickness  t  and  length  l  arranged  parallel  to  the  wall  (i.e.  at  zero  angle  of  attack),  offset  by  some 
height  h  <  do.  LEBU  configurations  with  more  than  one  parallel  plate,  sometimes  referred  to  as 
multiple  element  LEBUs,  are  shown  in  Figure  2.1b  and  Figure  2.1c.  Arrangements  of  vertically 
stacked  parallel  plates  of  identical  length  and  thickness,  later  referred  to  as  Multi  LEBU  devices, 
are  shown  in  Figure  2.1b,  along  with  the  notation  used  to  identify  the  height  of  each  plate.  For 
multi-LEBU  devices,  the  heights  of  each  plate,  starting  with  the  outermost  and  proceeding 
towards  the  wall,  are  denoted  as  h\,  and  hi.  LEBU  plates  arranged  in  series  in  the  streamwise 
direction  are  later  referred  to  as  Tandem  LEBU  devices,  as  shown  in  Figure  2.1c.  Each  individual 
plate  is  identical  in  length  /,  and  is  mounted  at  the  same  height  off  the  wall  h,  but  plates  are 
separated  by  some  gap  distance  5  in  the  streamwise  direction. 

Effect  of  LEBU  Devices  on  TBL  Characteristics 

The  characteristics  of  LEBU  modified  boundary  layers  have  been  well  documented  using 
an  array  of  diagnostic  methods  in  order  to  quantify  and  maximize  viscous  drag  reductions,  and  to 
investigate  the  mechanisms  that  were  responsible  for  LEBU  effects  (Anders,  1990,  and  the 
references  therein).  Smoke-wire  flow  visualization  records  (Corke,  1981;  Savill  &  Mumford, 
1988)  showed  a  significant  change  in  turbulence  structure  in  the  outer  part  of  the  LEBU  modified 
TBL.  In  particular,  the  flow  visualization  images  indicate  that  the  size  of  turbulent  structures  in 
the  outer  layer  are  substantially  reduced  (inspiring  the  LEBU  moniker),  and  showed  an  obvious 
reduction  in  intermittency  in  the  TBL  edge,  which  would  reduce  the  entrainment  of  irrotational, 
high  momentum  freestream  fluid  (Corke,  1981). 

Potential  for  Aero-Optic  Mitigation 

Prior  experimental  work  exploring  the  aero-optics  of  boundary  layers  has  shown  strong 
evidence  for  the  dominance  of  large-scale  (i.e.  outer-region)  turbulent  structures  as  a  source  for 
boundary  layer  induced  aberrations.  It  is  also  suspected  that  large  instantaneous  drop-outs  in  far- 
field  beam  intensity  caused  by  TBLs  are  related  to  large-scale  motions  (LSMs)  in  the  outer  part  of 
the  TBL  (Gordeyev,  Cress,  &  Jumper,  2013)  such  as  turbulent  bulges  and  hairpin  vortex  (Cress, 
2010).  Conveniently,  LEBU  devices  are  effective  at  reducing  streamwise  and  spanwise  integral 
length  scales,  boundary  layer  intermittency,  and  the  strength  of  bulges  and  sweeps  without 
introducing  a  strong  shear  layer  or  other  optically  un-desirable  flow  features,  making  them  an 
attractive  candidate  for  use  as  an  aero-optic  mitigation  application  for  TBLs. 

While  the  potential  for  using  LEBUs  for  aero-optic  mitigation  was  recognized  over  three 
decades  ago  by  Jumper  (1980),  and  again  several  years  later  by  Anders  (1985),  these  suggestions 
have  not  been  seriously  investigated.  It  is  important  to  recall  that  prior  studies  of  mechanisms  and 
optimization  of  LEBU  devices  have  been  focused  on  reducing  skin  friction,  which  is  related  to  the 
condition  of  the  TBL  locally  near  the  wall.  Also,  the  condition  of  having  a  small  amount  of  LEBU 
device  drag,  which  was  an  important  factor  for  LEBUs  drag  reduction  optimization,  is  not 
necessarily  a  restrictive  factor  for  aero-optic  mitigation.  In  real  directed  energy  applications,  it 
may  be  desirable  to  accept  a  small  additional  localized  drag  penalty  in  exchange  for  improved 
aero-optic  performance  over  a  finite  aperture. 

2.2  Experimental  Details 

Experimental  results  presented  in  this  report  were  performed  in  the  Transonic  Wind 
Tunnel  (TWT)  facility  located  in  Hessert  Laboratory  for  Aerospace  Research  at  the  University  of 
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Notre  Dame.  The  TWT  is  composed  of  an  inlet  contraction,  with  an  area  ratio  of  150:1  and 
screens  and  honeycombs  to  reduce  ffeestream  turbulence  intensities  (measured  w/  hot-wires  to  be 
on  the  order  of  2%).  Following  the  inlet  contraction  is  a  modular  smooth  wall  boundary  layer 
development  section  10.0  cm  x  9.9  cm  (or  about  4"x4"),  an  optical  measurement  section,  and  a 
diffuser  section  that  leads  to  the  vacuum  pump  plenum. 

In  the  present  study,  the  boundary  layer  development  length  to  the  optical  measurement 
location  varied  between  100  cm  and  170  cm.  In  the  optical  measurement  section,  the  upper  and 
lower  Plexiglas  walls  were  replaced  with  optical  quality  glass  in  order  to  ensure  good  optical 
access  for  aero-optic  wavefront  measurements.  Freestream  Mach  number  spanned  the  range  from 
approximately  0.2  to  0.5.  The  boundary  layer  thickness,  8,  was  found  to  be  1.9  cm,  and  the 
integral  displacement  and  momentum  thicknesses  were  2.6  mm  and  1.9  mm,  respectively.  At  170 
cm,  8  was  found  to  be  approximately  2.4  cm,  8  was  3.6  mm,  and  the  momentum  thickness  0  was 
2.75  mm.  The  shape  factors  for  both  of  these  cases  ranged  from  H  =  1.31  -  1.33,  which  is 
consistent  with  other  TBLs  at  these  Reynolds  numbers  (Nagib,  et  al.  2007). 

The  two-beam  Malley  probe  used  in  this  work  used  a  20  to  30  mW  continuous  wave  (CW) 
HeNe  (2  =  633  nm)  laser  source  that  was  spatially  filtered  to  generate  a  collimated  small  diameter 
beam  with  high,  near  uniform  intensity  across  the  beam  profde.  The  spatial  filter  was  also  used  to 
control  the  diameter  of  the  source  beam,  which  was  approximately  1  mm.  The  source  beam  was 
then  split  into  two  parallel  beams  of  equal  intensity  using  a  50/50  beam  splitter  plate  and  a 
steering  mirror.  The  spacing  between  these  beams,  A,  is  typically  on  the  order  of  5  -  10  mm.  This 
pair  of  parallel  beams  then  passes  straight  through  a  beam  cube  (losing  50%  of  its  intensity)  and 
through  the  measurement  region  (typically  a  wind  tunnel  test  section)  to  the  return  mirror.  The 
return  mirror  is  aligned  so  that  the  Malley  beams  are  reflected  back  to  the  beam  cube  along  the 
same  path.  If  the  deflection  angles  are  small  (which  they  are  for  aero-optic  measurements)  then 
this  is  effectively  an  analog  doubling  of  the  signal  intensity  for  the  Malley  probe. 


Figure  2.2:  Schematic  of  a  typical  two-beam  Malley  probe  sensor  in  a  double-pass,  double 
boundary  layer  (DBL)  wavefront  measurement  configuration. 

The  beam  cube  reflects  the  return  beams  onto  a  configuration  of  turning  mirrors  that 
isolate  each  of  the  two  parallel  Malley  beams.  Each  beam  is  then  passed  through  a  final  focusing 
lens,  which  in  this  case  has  a  focal  length  of  1  m,  onto  an  analog  position  sensing  device  (PSD) 
placed  in  the  focal  plane  of  the  final  focusing  lens.  The  PSD  measures  the  displacement  of  the 
beam  spot  location  on  its  surface.  Since  the  PSD  is  located  one  focal  length  away  from  the  final 
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focusing  lens,  a  ray  incident  on  the  lens  that  is  deflected  at  some  angle  0  will  can  be  calculated 
from  the  beam  displacement  on  the  PSD. 

2.3  Results  and  Discussion 

Three  different  LEBU  geometries  were  tested:  a  single  LEBU,  Figure  2.1a)  and  multiple  element 
LEBU  devices  called  Multi  LEBU,  Figure  2.1b)  and  Tandem  LEBU,  Figure  2.1c).  All  relevant 
LEBU  parameters  are  also  defined  in  Figure  2.1.  One  goal  of  this  investigation  is  to  verify  that  the 
LEBUs  used  in  this  work  modify  the  velocity  field  in  a  manner  that  is  consistent  with  the  results 
of  prior  studies.  The  present  work  also  study  seeks  to  study  the  relationship  between  the 
characteristics  of  the  LEBU-modified  TBLs  and  their  corresponding  effect  on  wavefront 
statistics.  It  also  presents  the  results  of  parametric  studies  of  single  and  multiple  element  LEBU 
devices  in  order  to  determine  the  optimal  LEBU  configuration  for  aero-optic  mitigation. 

2.3.1  Single  LEBU  Configuration 

A  full  record  of  the  LEBU  configurations  that  were  investigated  in  this  work  is  presented 
in  Table  2.1.  During  all  experimental  measurements,  ffee-stream  velocity  was  held  constant  for  all 
tests  at  approximately  M  =  0.4. 


Table  2.1.  Single  LEBU  Device  Configurations 


Length,  US 

0.8 

1.0 

1.6 

4.0 

5 

0.3 

♦ 

♦ 

♦ 

♦ 

< 

0.5 

♦ 

♦ 

♦ 

♦ 

X 

JDJD 

0.6 

♦ 

♦ 

❖ 

♦ 

X 

0.8 

♦ 

♦ 

♦ 

♦ 

Legend 

♦ 


Malley  probe  wavefront  measurements  single  modified  TBL 


Single  and  double  modified  TBL  cases  for  wavefront  measurements 


Velocity  Measurements 

To  investigate  the  effect  of  the  devices  on  the  turbulent  boundary  layer,  hot-wire  velocity 
measurements  were  obtained  at  three  streamwise  locations  (x/S  =  1.5,  3.5  and  5.4)  downstream  of 
the  single  LEBU  devices  using  the  hot-wire  anemometer.  From  these  data,  mean  and  fluctuating 
velocity  profiles  were  computed,  along  with  power  spectral  densities  of  fluctuating  velocity,  and 
compared  to  results  of  the  baseline  TBL  in  order  to  evaluate  the  effect  of  LEBU  devices  on  the 
boundary  layer.  These  data  were  also  compared  to  prior  LEBU  studies  in  order  to  confirm  that  the 
LEBU  devices  used  in  this  experiment  behave  similarly  to  what  previous  authors  have  reported. 

Mean  velocity  profiles  are  plotted  in  Figure  2.3  for  hot-wire  velocity  measurements 
obtained  at  three  streamwise  locations  (x/S  =  1.5,  3.5,  and  5.4)  downstream  of  the  LEBU  trailing 
edge  for  of  the  /  =  1.6<S,  h  =  0.6 S  LEBU  device.  The  wall-normal  coordinate  y  is  normalized  by 
the  incoming  TBL  thickness  measured  at  the  LEBU  location,  So  =  2.4  cm.  The  velocity  profile  for 
the  un-perturbed  boundary  layer  is  plotted  for  reference,  and  a  horizontal  dashed  line  marks  the 
LEBU  height.  Downstream  of  the  LEBU  devices,  a  well-defined  wake  deficit  is  observed 
imprinted  on  the  TBL  mean  velocity  profiles;  this  wake  region  is  centered  about  the  LEBU 
height,  and  as  x!5  increases  the  maximum  velocity  deficit  decreases  and  the  wake  spreads  in  the 


9 


spanwise  direction.  Outside  of  the  LEBU  wake  region,  there  is  no  significant  change  observed  in 
the  mean  velocity  profile. 


Figure  2.3.  Mean  velocity  profiles  of  the  LEBU-modified  TBL  at  several  locations  downstream  of 
a  Single  LEBU,  /  =  1.6 S,  h  =  0.6 3. 


M rms  (y)/ux 

Figure  2.4.  Root-mean-square  velocity  profiles  of  the  LEBU-modified  TBL  at  several  locations 
downstream  of  a  Single  LEBU,  l  =  \.6d,  h  =  0 .63. 

Profiles  of  the  root-mean-square  of  the  fluctuating  velocity  component  for  the  same  LEBU 
device,  which  are  plotted  in  Figure  2.4  along  with  the  baseline  case,  showed  a  significant 
reduction  in  urms  in  the  LEBU  modified  boundary  layer,  with  the  maximum  reduction  being 
located  at  the  LEBU  height.  Unlike  the  mean  velocity  profiles,  the  magnitude  of  the  reduction  did 
not  change  significantly  with  increasing  x,  but  rather  remained  constant  around  urms  ~  0.04 U®. 
Velocity  fluctuations  in  the  region  of  the  boundary  layer  below  the  LEBU  height  also  showed  a 
reduction  in  magnitude  with  respect  to  the  baseline,  with  the  reductions  extending  further  towards 
the  wall  with  increased  distance  downstream  of  the  LEBU  trailing  edge.  Above  the  LEBU  device, 
small  peaks  were  found  in  the  urms  profiles  at  all  streamwise  locations  that  exceeded  the  baseline 
Urms  locally.  The  magnitude  of  these  spikes  in  RMS  decreased  further  downstream  from  the  LEBU 
and  the  peak  location  migrated  away  from  the  wall  slightly. 
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Upon  inspection,  it  was  found  that  these  velocity  data  are  consistent  with  typical  mean 
velocity  profiles  obtained  in  previous  studies  (Corke,  et  al.,  1979;  Hefner,  et  al.,  1979;  Lemay,  et 
al.,  1990),  which  showed  a  similar  velocity  deficit  in  the  region  just  downstream  of  LEBU 
devices. 

Wake-Deficit  Velocity,  Profiles 

In  the  mean  velocity  profiles  shown,  the  maximum  velocity  deficit  in  the  wake  region  is 
shown  to  diminish  as  downstream  distance  increases,  indicating  that  the  modified  boundary  layer 
is  recovering  from  the  disturbance  caused  by  the  LEBU  device,  at  least  in  the  mean  velocity 
profiles.  The  spanwise  extent  of  the  wake  imprint  spreads  in  the  wall-normal  direction  with 
increasing  distance  downstream.  Both  of  these  characteristics  are  consistent  with  the  behavior  of 
planar  self-similar  wakes  in  constant  mean  flow  like  those  investigated  by  Moser,  et  al.  (1998). 
They  found  that  for  planar  wakes,  the  wake  velocity  deficit,  Uo,  and  the  wake  half-width,  b, 
evolved  in  the  streamwise  direction  as  b  ~  xm  and  Uo  ~  x’1/2.  The  non-dimensional  growth  rate  of 
the  wake,  defined  as  a  =  0_1  d(b2  )/dx ,  where  0  is  the  wake  momentum  thickness,  was  found  to 
be  0.29  -  0.41  for  canonical  self-similar  plane  wakes. 

To  determine  if  the  LEBU  wake  evolves  in  a  self-similar  manner  within  the  TBL,  a 
modified  form  of  the  wake  deficit  velocity  profile,  U Deficit  (>’),  is  computed  by  subtracting  the 
LEBU-modified  boundary  layer  profiles  from  the  baseline  velocity  profiles  (Moser,  et  al.,  1998): 

U Deficit  GO”  ^ LEBU  GO”  U Baseline  GO  (2.1) 

The  maximum  wake  deficit  velocity,  Uo,  is  the  maximum  of  the  absolute  value  of  the  wake  deficit 

profile,  max|^£>ey!Cft(f)|J  Note  that  where  the  baseline  velocity  is  greater  than  the  mean 
velocity  of  the  LEBU  profile,  UDefidt  will  be  negative.  The  wake  half-width,  b,  is  defined  as  the 
distance  between  the  two  points  where  -  U  Deficit  (y)  =  0.5 U0  (Moser,  et  al.,  1998). 


Figure  2.5.  Mean  wake  velocity  deficit  profiles  downstream  of  the  /  =  1.6 5,  h  =  0.6b  single  LEBU 
device. 
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Mean  velocity  deficit  profiles  for  the  LEBU  modified  boundary  layer  data  were  computed, 
and  the  results  are  shown  in  Figure  2.5  for  the  /  =  1 .63,  h  =  0.63  LEBU  device,  along  with  the 
mean  deficit  profile  from  Moser,  et  al.  (1998)  for  canonical  plane  wakes.  The  LEBU  wake  deficit 
profiles  presented  exhibit  self-similarity  that  is  consistent  with  canonical  plane  wake  data  on  the 
outer  part  of  the  boundary  layer  (i.e.  (y-  h)  lb>  0)  for  all  streamwise  locations.  For  the  inner  part 
of  the  wake  (i.e.  (y  -  h)  /b<  0),  the  wake  profiles  deviate  slightly  from  the  self-similarity  scaling, 
with  the  deviation  becoming  larger  for  locations  further  downstream.  From  these  data,  the  non- 
dimensional  growth  rate  was  found  to  be  approximately  0.15  for  the  LEBU  modified  TBL  wake 
profiles.  This  was  much  less  than  the  range  of  growth  rates,  0.29  -  0.41,  found  for  canonical  self- 
similar  plane  wakes  (Moser  et  al,  1998). 

The  deficit  profile  for  the  RMS  velocity  fluctuations  can  be  computed  in  a  similar  fashion 
to  the  mean  velocity  in  equation  (2.1): 

^ rms ,  Deficit  0^ )  —  ^  rms,  LEBU  0^ )  ^ rms ,  Baseline  (y)  (2.2) 

RMS  velocity  deficit  profiles  computed  via  equation  (2.2),  however,  were  not  found  to  scale  well 
with  self-similar  planar  wake  data,  due  to  the  large  reductions  in  urms  below  the  LEBU  device. 
The  location  of  the  ‘spike’  in  unns  that  occurs  on  the  upper  side  of  the  LEBU  in  Figure  2.4  is 
found  to  be  constant  at  about  (y~h)  Id-  0.6.  This  is  an  indication  that  the  local  increase  in  urms 
above  the  LEBU  is  closely  related  to  the  wake  deficit  profile  in  the  modified  TBL. 


(y  ~  h)/b 

Figure  2.6.  Wake  root-mean-square  velocity  profiles  downstream  of  the  /  =  1.63,  h  =  0 .63  single 
LEBU  device.  Legend  is  the  same  as  Figure  2.5. 

Corke  (1981)  recognized  that  local  changes  in  the  mean  velocity  profile  could  have  a 
significant  impact  on  velocity  fluctuations  in  the  LEBU-modified  boundary  layer,  since  the 

dominant  turbulence  production  term,  u'v\dU  /  dy) ,  is  directly  proportional  to  the  local  mean 
velocity  gradient.  Since  the  wake  deficit  mean  profile  induced  on  the  LEBU  within  a  sheared 
mean  flow,  the  local  velocity  gradient  below  the  LEBU  device  is  decreased  and  the  velocity 
gradient  above  the  LEBU  is  increased.  This  is  illustrated  in  Figure  2.7,  where  dU/dy  and  the 
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fluctuating  velocity  component  in  the  LEBU-modified  boundary  layer  are  plotted  for  easy 
comparison.  Here,  dU/dy  was  computed  from  experimental  data  using  a  first  order  central 
difference  method.  It  is  evident  from  Figure  2.7  that  there  is  some  correlation  between  the 
velocity  gradient  and  urms,  especially  in  the  neighborhood  of  the  LEBU  height,  however  changes 
in  dU/dy  for  the  LEBU-modified  boundary  layer  does  not  fully  account  for  reductions  in  velocity 
fluctuations  below  the  LEBU  device.  Rather  it  is  likely  that  both  the  changes  in  velocity  gradient 
and  the  limitation  of  vertical  velocity  fluctuations  (i.e.  the  plate  effect),  that  suppressed  urms  over 
such  a  wide  range  of  the  area  below  the  LEBU  device.  Additional  support  for  this  notion  was 
shown  in  Lemay,  et  al.  (1990),  where  measurements  showed  suppressed  production  for  y  <  h  in 
the  LEBU-modified  boundary  layer  up  through  the  location  where  the  mean  profile  (and  therefore 
the  velocity  gradient)  is  nearly  recovered  to  baseline  state. 

a)  b) 


7/2  _ (».Baseline\2  Q 

“'rms  \arms  )  ^  10'1 


oo 


Figure  2.7.  a)  Normalized  mean  velocity  gradient  and  b)  baseline-normalized  urms  profiles  for 
LEBU  modified  boundary  layer  (1/3  =  1.6,  hi 5  =  0.6).  Legend  is  the  same  as 
Figure  2.4. 


Deflection  Angle  Spectra 

Amplitude  spectra  of  Malley  probe  deflection  angle  measurements  were  computed  for 
measurements  downstream  of  the  single  LEBU  devices  at  a  number  of  locations  up  to  lOd 
downstream  of  the  device  trailing  edge.  The  results  are  presented  for  several  streamwise  locations 
in  Figure  2.8  where  frequency  is  non-dimensionalized  as  Strouhal  number  based  on  freestream 
velocity,  Vx,  and  the  canonical  TBL  thickness  at  the  LEBU  mounting  location,  3.  Note  that  for  St„ 
<0.2,  deflection  angle  spectra  are  contaminated  by  mechanical  vibration  from  the  wind  tunnel 
motors.  Also,  in  the  high  frequency  range  of  the  spectra  the  sharp  peaks  near  Sf>  =  6,  7,  and  10 
were  found  to  be  the  result  of  electronic  interference  from  sources  within  the  laboratory 
environment. 
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St* 


St.* 

Figure  2.8.  Deflection  angle  spectra  from  Malley  probe  measurements  of  the  modified  TBL  for 
the  l  =  1.6(5,  h  =  0.6(5  single  LEBU  device  at  streamwise  locations  a)  from  x  =  0.8(5  to  x  =  4.1(5  and 
b)  from  x  =  4. 1(5  to  x  =  9.5(5. 

Figure  2.8  presents  deflection  angle  spectra  obtained  from  Malley  probe  measurements 
downstream  of  the  /  =  1.6(5,  h  =  0.6(5  single  LEBU  device  at  several  streamwise  locations.  Just 
downstream  of  the  device  (x  ~  Id),  the  spectra  presented  in  Figure  2.8a  showed  significant 
reductions  over  the  range  of  low  frequencies,  0.2  <  Stg  <  4.  Spectra  at  streamwise  locations  x<  4(5 
did  not  show  any  significant  increase  over  the  baseline  deflection  angle  spectrum;  by  x  =  2.4(5  the 
peak  of  the  LEBU-modified  TBL  spectrum,  which  is  around  SL  =  2,  increased  slightly.  The  peak 
amplitude  then  remained  relatively  constant  through  approximately  x  =  4d.  Portions  of  the  high- 
frequency  range  of  the  spectra  (SL  >  4)  and  the  low-frequency  end  of  the  spectra  (St,,  <  0.5)  were 
found  to  decrease  with  increasing  x  between  3-4(5  downstream.  Overall,  for  this  LEBU 
configuration,  the  reductions  in  the  deflection  angle  spectra  were  found  to  be  the  largest  compared 
to  other  tested  single  LEBU  devices. 

Deflection  angle  spectra  collected  further  downstream  between  x/d  =  4  -  10  for  the  h  = 
0.6d  single  LEBU  are  shown  in  Figure  2.8b.  In  general,  spectra  from  these  distances  downstream 
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of  the  LEBU  trailing  edge  were  found  to  result  in  broadband  reductions  that  were  notably  larger 
than  what  was  found  for  other  single  LEBU  devices.  The  peak  amplitude  of  the  deflection  angle 
was  shown  to  be  increasing  for  measurements  further  downstream  of  the  LEBU,  and  the  peak 
location  was  found  to  be  shifting  slightly  back  towards  St<5  =  1,  the  peak  location  for  the  un¬ 
perturbed  TBL  wavefront  spectra.  As  x/S  approached  10,  there  were  notable  increases  in  the  low- 
frequency  range  of  the  spectra  that  were  caused  by  the  shift  in  peak  frequency  and  increase  in 
peak  amplitude  of  the  LEBU-modified  TBL. 

In  summary,  inspection  of  deflection  angle  spectra  for  the  single  LEBUs  presented  in 
Figure  2.8  and  other  tested  ones  demonstrated  that  the  devices  yielded  immediate  and  significant 
reductions  in  the  low-frequency  range  of  the  spectra.  Note  that  assuming  frozen  flow  (which  has 
been  shown  to  be  a  good  assumption  for  Malley  probe  wavefront  measurements  in  the  TBL),  the 
range  of  frequencies  over  which  reductions  were  observed  corresponds  to  turbulent  structures  that 
are  up  to  53  in  length.  For  Malley  probe  spectra  for  the  LEBU  modified  TBLs,  the  peak  location 
was  found  to  initially  shift  to  higher  Strouhal  numbers  (~  2),  indicating  that  the  characteristic 
length  scales  in  the  modified  TBL  are  reduced  via  the  LEBU  devices  and  their  effect  on  large- 
scale  motions. 

The  behavior  of  the  high-frequency  portion  of  the  deflection  angle  spectra  gave  some 
evidence  that  the  energy  formerly  contained  in  the  large  scale  turbulent  structures  was 
redistributed  by  the  LEBU  wake  into  smaller  scale  turbulent  structures.  The  turbulent  kinetic 
energy  (TKE)  contained  in  these  small  scale  structures  was  quickly  dissipated  through  the  energy 
cascade.  After  the  excess  TKE  is  shed  through  the  mechanism  of  dissipation,  the  high-frequency 
range  of  spectra  is  also  reduced.  These  observations  are  consistent  with  measurements  of  the  TKE 
balance  in  a  LEBU  modified  TBL  performed  by  Lemay,  et  al.  (1990),  which  showed  both  a 
significant  increase  in  the  production  term  and  an  increase  in  the  dissipation  term  in  the  LEBU 
wake  in  the  neighborhood  of  the  y  =  h  just  downstream  of  the  LEBU  device  (approx.  x/S  ~  2.5). 
Within  10d  downstream  of  the  LEBU  trailing  edge,  however,  Lemay,  et  al.  (1990)  found  that  the 
dissipation  term  is  nearly  recovered  to  the  baseline  canonical  TBL  levels.  This  is  consistent  with 
the  downstream  distances  found  in  the  present  study  at  which  wavefront  spectra  were  found  to 
show  reductions  in  the  high-frequency  end  of  the  spectra. 

To  better  visualize  the  evolution  of  deflection  angle  spectrum  in  the  streamwise  position 
as  compared  to  the  baseline  deflection  angle  spectrum,  the  ratio  of  deflection  angle  spectrum  for 
LEBU-modified  TBLs  to  baseline  TBL  spectrum  can  be  defined  as 


Ce{x/5,Sts ) 


^LF.BU  (X/3  5  Stg  ) 

^Baseline  (X/^  ’  ) 


(2.3) 


where  Co  (x/d,  Sts)  will  be  equal  to  zero  where  the  LEBU-modified  spectrum  is  equal  to  the 
baseline  spectrum,  positive  where  the  LEBU-modified  spectrum  is  increased  compared  to  the 
baseline,  and  negative  where  the  LEBU-modified  spectrum  is  less  than  the  baseline  spectrum. 
Contour  plots  of  Co  (xJS,  Sts)  are  presented  in  Figure  2.9for  both  single  LEBU  devices.  For  both 
cases,  Figure  2.9  shows  a  significant  reduction  in  large-scale  structures  below  Stg  =  3  that 
diminishes  as  distance  from  the  LEBU  trailing  edge  increases.  Additionally,  a  decrease  in  high- 
frequency  content  of  the  deflection  angle  spectrum  is  also  shown  as  the  downstream  distance 
increases,  which  is  consistent  with  the  results  previously  shown  in  Figure  2.8. 
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Figure  2.9.  Baseline-normalized  spectrum  surfaces  Co  (x/d,  Sts )  for  single  LEBU  devices  of  height 
a)  hid  =  0.8  and  b)  hid  =  0.6.Note  that  *  denotes  the  peak  location  of  Co  at  each  value  of  x. 


One  feature  that  is  readily  apparent  in  the  LEBU-modified  deflection  angle  spectrum 
plotted  in  Figure  2.9  is  the  evolution  of  the  peak  of  Co  (x/d,  Sts )  with  the  downstream  distance 
from  the  LEBU  device,  which  likely  corresponds  to  the  development  of  the  LEBU  device  wake. 
The  shift  of  the  peak  in  Co  to  the  lower  frequencies  occurs  as  the  downstream  distance  increases, 
eventually  approaching  a  value  of  Sts  =  L  Several  authors  (Corke,  1981;  Lemay,  et  al.  1990,  and 
others)  found  that  just  downstream  of  the  LEBU  there  is  a  momentum-deficit  wake  that  begins  at 
the  LEBU  trailing  edge,  and  grows  with  the  increasing  downstream  distance,  x,  to  a  point  where  it 
is  on  the  order  of  the  TBL  thickness  d  and  interacts  with  the  near-wall  layer.  If  it  is  assumed  that 
the  peak  in  Co  corresponds  to  a  coherent  structure  of  some  size,  y,  then  the  relationship  between 
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the  structure  size  y  and  the  frequency  fpeak  at  which  the  maximum  in  Co  occurs  can  be  expressed 
via  the  frozen  flow  assumption, 


r  = 


u„ 


f, 


Peak 


8 

Stpseak 


(2.4) 


where  St[eak  is  the  Strouhal  number  that  corresponds  to  fPeak-  If  it  is  assumed  that  the  shift  of 
St/**  toward  lower  frequencies  corresponds  to  the  growth  in  the  size  of  the  LEBU  wake,  then  the 
half-width,  b,  of  the  LEBU  wake  can  be  estimated  as  a  function  of  x  from  the  peak  of  the 
normalized  spectrum; 


b{x) 


1  S 

2  StPeak 


(2.5) 


The  result  of  this  estimation  for  the  single  LEBU,  h  =  0.63  device  is  presented  in  Figure  2.10 
along  with  measurements  of  the  wake  half-width  computed  from  hot-wire  velocity  profiles.  The 
comparison  shows  that  the  estimate  of  the  wake  half-width,  b,  from  the  normalized  spectrum 
peaks  is  in  good  agreement  with  direct  measurements  obtained  previously  for  the  same  case,  and 
verifies  that  the  evolution  of  the  peak  of  the  normalized  spectrum  Co  is  a  result  of  the  growth  of 
the  LEBU  wake  within  the  modified  TBL. 
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Figure  2.10.  Comparison  of  wake  half-width  b  estimated  from  Co  (x/S,  Sts )  and  from  hot-wire 
velocity  measurements  downstream  of  the  /  =  1.6<5,  h  =  0 .63  from  section  6.2.1. 
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2.3.2  Streamwise  Wavefront  Correlations 

The  streamwise  wavefront  correlation  functions  of  Malley  probe  wavefront  data  were 
computed  by  taking  the  inverse  Fourier  transform  (Fl)  of  the  single-boundary  layer  wavefront 
frequency  spectra  to  obtain  the  autocorrelation  in  time  at  each  streamwise  measurement  location: 

Rw{r)  =  F-l\w{f)  |2}.  (2'6) 

The  frozen-flow  approximation,  Ax  =  Ucz,  was  then  applied  to  estimate  the  streamwise  wavefront 
correlation  function,  Rw  (Ax).  The  correlation  coefficient  was  then  defined  by  normalizing  Rw  by 
the  peak  value: 
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(2.7) 


Pw  (A^) 


(a^) 
Rj 0)  ' 


To  calculate  the  integral  correlation  length,  the  correlation  coefficient  pw  (x,  Ax)  was  integrated 
along  the  Ax  coordinate  to  the  location  of  the  first  zero  crossing,  Axzero: 

Ax 

A  W(x/S)  =  j  pw(x,  Ax)d(Ax) 
o 

This  choice  was  made  because  the  negative  correlation  values  past  the  first  zero  crossing  are  not 
large  in  magnitude;  it  has  also  been  shown  previously  that  this  region  of  the  correlation  function 
is  highly  sensitive  to  finite  aperture  effects  (Smith,  et  al.,  2012;  Gordeyev,  et  al.  2014). 

The  integral  correlation  length  was  computed  for  the  baseline  TBL  wavefront  data  and  for 
both  single  LEBIJ  cases.  The  results  for  both  single  LEBIJ  devices  are  shown  in  Figure  2.14, 
normalized  by  the  baseline  correlation  lengths,  in  order  to  determine  the  relative  reduction  in 
streamwise  correlation  length  for  the  modified  TBLs.  For  both  devices  these  results  show  definite 
reductions  of  15  -  34%  in  correlation  length  for  the  LEBU  modified  TBLs,  depending  on  the 
LEBU  configuration  and  streamwise  location.  LEBU  with  h  =  0.6<5  gives  the  largest  reduction  in 
streamwise  correlation  length  over  the  whole  of  the  optical  measurement  section,  with  a 
maximum  reduction  of  approximately  34%  at  x  =  4  -  5d. 


Figure  2.1 1.  Integral  correlation  length  computed  from  equation  (2.8)  as  a  function  of  streamwise 
location. 


2.3.3  Single  LEBU  Parametric  Study 

In  order  to  identify  the  optimal  single  LEBU  device  parameters  that  yield  the  largest 
sustained  reductions  in  OPD ^  across  the  optical  measurement  section  in  the  present  study,  a 
parametric  investigation  of  LEBU  chord  length,  l,  and  height,  h,  was  performed.  For  this 
investigation,  the  LEBU  length  was  varied  between  /  =  0.8 3  and  4 <5,  and  the  height  was  varied 
between  h  =  0.3d  and  either  0.8<5  or  l.Od. 


Effect  of  Device  Length  and  Height 

To  investigate  the  effect  of  LEBU  device  length  on  OPDrms  in  the  manipulated  turbulent 
boundary  layer,  Figure  2.12  presents  the  results  of  wavefront  measurements  obtained  at  fixed 
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height  h  =  0.63  for  several  different  chord  lengths.  This  particular  height  was  chosen  because  it 
gave  the  best  results  for  the  preliminary  single  LEBU  investigation  presented  in  the  previous 
chapter.  Note  that  these  data  are  normalized  by  the  baseline  OPDrms,  which  was  measured  for  the 
un-modified  TBL. 


x/60 

Figure  2.12.  Baseline-normalized  OPDrms  measured  downstream  of  single  LEBU  devices  with 
fixed  height  (h  =  0.63)  for  different  chord  lengths,  /. 

It  is  evident  from  the  data  presented  in  Figure  2.12  that  device  chord  length,  /,  had  a 
significant  impact  on  the  level  of  aero-optic  mitigation  that  was  achieved.  In  general,  longer 
devices  yield  larger  reductions  than  shorter  devices;  the  longest  LEBU  device  tested  was  4d  in 
length,  and  gave  reductions  in  OPDrms  on  the  order  of  40  -  45%  over  a  significant  portion  of  the 
measurement  section.  Based  on  these  results,  it  is  possible  that  longer  LEBUs  could  be  even  more 
effective  at  suppressing  aero-optic  aberrations  in  the  TBL.  This  result  is  consistent  with 
experimental  findings  from  Plesniak  (1984)  and  Savill  &  Mumford  (1988)  for  C/ reductions  in 
LEBU  manipulated  TBLs,  who  tested  LEBUs  up  to  5b  in  length  and  found  no  global  optimum 
value  for  1/3. 


Figure  2.13.  Baseline-normalized  OPDrms  measured  downstream  of  a  long  single  LEBU  (/  =  43) 
at  different  plate  heights,  h. 
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Baseline-normalized  OPDrms  obtained  for  the  long,  single  LEBIJ,  /  =  4S  device  at  several 
different  heights  are  presented  in  Figure  2.13  in  order  to  demonstrate  the  effect  of  LEBU  height 
on  the  aero-optic  characteristics  of  the  manipulated  TBL.  For  the  case  where  the  LEBU  was  the 
closest  to  the  wall,  h  =  0.3 3,  there  was  about  a  15%  reduction  in  OPDrms  just  downstream  of  the 
device  trailing  edge.  Between  x  =  3  -  4 3,  OPDrms  reached  a  minimum  value  between  0.7  -  0.75  of 
the  baseline  value,  OPD^EUNE .  Downstream  of  this  location,  OPDnns  increased  to  about 

0.85  OPD^eune  by  x  =  63,  and  to  nearly  0.9  0PD^EUNE  as  x  approached  1 0(5.  For  the  LEBU 

placed  at  the  next-to-lowest  height  in  this  study,  h  =  0.5(5,  OPDrms  was  initially  reduced  by  greater 
than  20%,  and  by  2 3  downstream  of  the  device  trailing  edge  it  rapidly  dropped  to  below 
0.6  OPD^EELINE  (i.e.  >40%  reduction).  OPD rms  remained  below  0.55  -  0.6  OPD^sSELINE  for  the 

region  x~  2  -  53.  Between  x  =  5  -  13,  however,  it  increased  to  0.7  OPD^sSELINE  and  remained 
between  0.7  -  0.75  OPD^sSELINE  until  the  end  of  the  measurement  section  around  10(5. 

For  the  case  in  which  the  LEBU  was  placed  slightly  higher  up  in  the  boundary  layer  at  h  = 
0.6 3,  OPDrms  reduced  to  approximately  0.56  OPD^EELINE  by  33  downstream,  and  OPDrms  remained 

between  0.55  -  0.6  OPD^sSELINE  (i.e.  a  40  -  45%  reduction  in  aberrations)  up  to  x~  8.5 3.  Beyond 

this  point,  OPDrms  increased  from  0.6  to  almost  0.7  as  x  approached  10(5,  indicating  that  the  aero- 
optic  environment  in  the  manipulated  TBL  might  have  been  beginning  to  recover  back  towards 
the  canonical  state.  The  final  tested  LEBU  height,  h  =  0.8(5,  showed  similar  reductions  in  OPDrms 
greater  than  30%  were  sustained  for  the  remainder  of  the  measurement  section,  with  a  local 
minimum  of  0.63  OPDE^sSEL,NE  around  x  =  13. 

From  these  data,  it  appears  that  near  the  wall  the  manipulated  flow  produces  only 
moderate,  short  lived  reductions  in  OPDrms,  while  devices  placed  near  the  center  of  the  TBL  (h  = 
0.5  -  0.6(5)  yield  the  largest  reductions  in  TBL  induced  wavefront  aberrations.  The  LEBU  placed 
slightly  above  the  middle  of  the  TBL  at  h  =  0.63  yielded  the  largest  sustained  streamwise  length 
of  reduction  in  OPDrms.  In  this  case  reductions  in  OPDrms  were  on  the  order  of  45%  over  nearly 
63,  and  were  greater  than  30%  for  about  93.  For  the  device  closest  to  the  TBL  edge  ( h  =  0.8(5), 
OPDrms  was  reduced  by  more  than  30%  for  about  8(5  in  the  streamwise  direction,  but  the 
maximum  level  of  reductions  was  less  than  that  observed  for  the  h  =  0.6(5  case. 

Optimal  SingleLERU Configuration 

The  results  of  the  parametric  study  of  LEBU  length  and  height  and  their  effect  on  aero- 
optic  aberrations  in  the  TBL  were  summarized  concisely  by  analyzing  a  few  key  criteria  for 
evaluating  the  performance  of  each  configuration: 

•  the  minimum  value  of  OPD^sBU  /  OPDBESEUNE  obtained  for  each  configuration,  and 

•  the  continuous  length  L  over  which  OPDrms  is  reduced  below  some  threshold  value. 

The  threshold  values  of  20%  and  30%  reductions  in  OPDrms  were  chosen  based  on  the 
observations  from  experiments,  thus  the  lengths  for  each  are  denoted  L  (2o%)  and  L  (30%), 
respectively. 

The  first  measure  of  LEBU  effectiveness  is  plotted  in  Figure  2.14  as  a  function  of  LEBU 
device  height  for  devices  of  different  length.  The  results  clearly  show  that  for  all  chord  lengths, 
the  maximum  reduction  in  OPDrms  occurs  where  the  LEBUs  are  placed  around  the  center  of  the 
TBL  near  h  =  0.5  -  0.63.  It  is  also  evident  from  this  data  that  the  magnitude  of  reductions  in 
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OP Drms  monotonically  increases  with  device  chord  length  regardless  of  LEBU  height,  with  the 
exception  of  US  =  0.8  and  1.0  at  h  =  0.8(5.  This  behavior  is  consistent  the  findings  of  parametric 
studies  from  previous  authors  (Plesniak,  1984;  Savill  &  Mumford,  1988)  investigating  local  C/ 
reductions  in  LEBU  manipulated  TBLs,  who  tested  LEBUs  up  to  5S  in  length  and  found  no 
global  optimum  value  for  l/S  with  regard  to  the  maximum  reduction  in  C/.  Recall  that  the  best 
value  often  cited  for  LEBU  drag-reduction  devices,  /  =  0.8  -  1.6(5  (Anders,  1990),  was  found  from 
optimizations  seeking  the  minimum  net  drag;  in  these  works  long  chord  lengths  were  found  to  add 
too  much  device  drag  to  be  effective. 


h/S 

Figure  2.14.  Minimum  values  of  OPD^u  /  OPD^sSELINE  as  a  function  of  device  height  for  single 
LEBU  devices  of  different  chord  lengths. 

Figure  2.15  presents  the  streamwise  lengths  L  (20%)  and  L  (30%),  over  which  OPDrms  is 
reduced  by  20%  and  30%  in  the  modified  TBL.  It  should  be  noted  here  that  streamwise  extent  of 
the  study  is  limited  to  10(5,  so  the  maximum  value  of  any  L  obtained  from  this  experiment  is  10(5. 
It  is  possible  that  for  some  single  LEBU  configurations  where  L  was  found  to  be  10<5  that  in 
reality  L  is  greater  than  1  Or)'.  The  present  study  refrains  from  speculation  as  to  which  cases  this 
may  or  may  not  occur  in,  but  an  extension  of  the  optical  measurement  section  beyond  10(5  in 
future  work  can  be  used  to  investigate  this  further  if  it  is  deemed  necessary.  In  Figure  2.15a  it  was 
found  that  the  continuous  length  over  which  OPDrms  was  reduced  by  20%  or  more  is  maximized 
for  h  =  0.6d  for  all  LEBU  devices.  LEBU  length  was  also  shown  have  a  notable  effect,  with  the 
value  of  L  (2o%)  increasing  with  US. 

The  results  for  L  (30%)  are  plotted  in  Figure  2.15b  as  a  function  of  LEBU  height,  for  devices 
of  different  length.  Note  that  increasing  the  threshold  to  reductions  greater  than  or  equal  to  30%  in 
OPDrms  only  preserves  cases  in  which  /  =  1.6-4 S  and  h  =  0.5  -  0.8(5.  Again  for  these  data,  h  = 
0.6<5  is  the  optimal  height  for  producing  long  regions  of  aero-optic  mitigation.  The  l  =  4S,h  =  0.6(5 
Long  LEBU  device  gives  the  largest  region  of  reduction  over  30%  L  (30%)  ~  9S.  Note  that  this 
configuration  also  provides  the  largest  local  reduction  in  OPDrms ,  as  shown  in  Figure  2.14. 
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Figure  2.15.  The  streamwise  distance  over  which  OPDrms  was  reduced  by  a)  more  than  20%  and 
b)  more  than  30%  with  respect  to  the  value  measured  for  the  un-manipulated  TBL. 


Table  2.2.  Multiple  Element  LEBU  Device  Configurations 

Multiple  Element  LEBU  Devices  (US  =  1.6) 


In/d 

h2/d 

s/d 

Measurement 

Type(s) 

Multi-LEBU 

0.9 

0.6 

- 

♦ 

0.8 

0.5 

- 

♦ 

0.6 

0.3 

- 

♦ 

Tandem  LEBU 

0.6 

- 

4.0 

♦ 

0.6 

- 

8.0 

❖ 

Legend 


♦ 


Malley  probe  wavefront  measurements  single  modified  TBL 


Single  and  double  modified  TBL  cases  for  wavefront  measurements 
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2.3.4  Multiple-Element  LEBU  Devices 

While  single  LEBUs  have  been  shown  to  give  sustained  reductions  in  OPDrms  on  the  order 
of  30%,  the  effectiveness  of  multiple-element  LEBU  devices,  including  multi-LEBU  (vertically 
stacked,  see  Figure  2.1b)  and  tandem  LEBU  (horizontally  spaced,  see  Figure  2.1c)  devices, 
should  be  assessed  to  investigate  whether  or  not  the  additional  elements  result  in  any  further 
reduction  in  levels  of  OPDrms. 

Multi-LEBU  Devices 

Overall  levels  of  OPDrms  were  computed  and  normalized  by  the  baseline  OPDrms,  these 
results  are  presented  in  Figure  2.16.  For  the  h\/6  =  0.9, h^d  =  0.6  and  h\!6  =  0.%,hiJ6  =  0.5  multi 
LEBU  configurations,  aero-optic  distortions  were  reduced  by  about  20%  immediately 
downstream  of  the  devices,  which  compares  well  with  the  reductions  observed  for  the  single 
LEBU  device  with  l/S  =  h/6  =  0.6. Further  downstream  of  the  LEBU  device,  the  levels  of  OPDrms 
for  both  of  these  multi  LEBU  devices  continue  to  decrease  with  respect  to  the  baseline 
measurements,  reaching  a  maximum  reduction  of  approximately  30%  around  the  streamwise 
position  of  16.  Compared  to  the  single  LEBU  result,  this  location  of  the  maximum  reduction  is 
approximately  26  further  downstream,  and  is  approximately  the  same  strength  as  the  reductions 
from  the  single  LEBU  device.  Beyond  the  location  of  the  maximum  reduction  for  both  of  these 
multi  LEBUs,  the  levels  of  OPDrms  for  h\/6  =  0.8,/72/b  =  0.5appear  to  begin  recovery  towards  the 
un-modified  boundary  layer  state,  while  OPDrms  for  h\/6  =  0.9, 62 16  =  0.6  is  shown  to  be  leveling 
out  around  70%  of  the  baseline  value.  These  results  appear  to  indicate  that  the  effect  for  multi 
LEBUs  do  not  seem  to  be  additive,  as  the  results  for  the  h\/6  =  0.9Jii/6  =  0.6  multi  LEBU  do  not 
show  a  simple  increase  in  the  reductions  detected  for  the  h!6  =  0.6single  LEBU  device. 


Figure  2.16.  Streamwise  development  of  OPDrms  for  multi  LEBU  devices  of  different  heights. 

For  the  h\/6  =  0.6, 6  =  0.3  multi  LEBU  configuration  OPDrms  was  reduced  by  about  30% 
immediately  downstream  of  the  device  trailing  edge,  and  OPDrms  continued  to  decrease,  resulted 
in  a  large,  relatively  flat  range  of  reduction  between  25  and  66  where  levels  are  about  40%  less 
than  those  of  the  baseline  OPDrms.  Beyond  66  downstream  of  this  multi  LEBU  device,  levels  of 
OPDrms  begin  to  recover  toward  the  un-modified  TBL  levels,  reaching  about  80  %  of  the  baseline 
value  by  lOd. 
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Figure  2.17.  Baseline-normalized  spectrum  surfaces  Ce  (x/d,  Sts)  for  two  Multi-LEBU  devices;  a) 
h\  =  0.9c),  hi  =  0.63,  and  b)  h\  =  0.63,  hi  =  0.33. 


Contour  plots  of  Ce  (x/d,  Sts),  which  was  defined  previously  in  equation  (2.3),  are 
presented  in  Figure  2.17  for  two  multi  LEBIJ  device  configurations.  For  the  h\/§  =  0.9, hit 3  =  0.6 
device,  shown  in  Figure  2.17a,  in  the  region  immediately  downstream  of  the  LEBU  (x  =  1  -43) 
there  is  a  significant  reduction  in  the  deflection  angle  spectrum  over  the  baseline  values  below  Sts 
~  2,  and  an  increase  above  the  baseline  values  in  the  high-frequency,  Sts>  3,  end  of  the  spectrum. 
This  increase  over  the  baseline  value  is  reduced  with  increasing  downstream  distance,  until  after  x 
~  43,  where  the  high-frequency  end  of  the  spectrum  (Sts>  2)  is  reduced  below  the  baseline  values. 
For  the  low-frequency  portion  of  the  deflection  angle  amplitude  spectrum  ( Sts<  2),  the  strength  of 
the  reduction  of  the  multi-LEBU  modified  spectrum  increases,  especially  around  a  value  of  Sts  = 
0.4. In  addition,  as  the  measurement  location  increases  from  4  to  lOd,  the  region  of  low-frequency 
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suppression  achieved  by  this  particular  LEBU  device  is  reduced  slightly  in  Strouhal  number  space 
as  the  peak  of  the  LEBU-modified  deflection  angle  spectrum  approaches  the  baseline  level  and 
baseline  location  around  Sts~  1 . 

The  normalized  spectrum  surface,  shown  in  Figure  2.17b,  for  the  h\/8  =  0.6,hi/8  =  0.3 
multi  LEBU  device  bears  some  qualitative  similarities  to  the  spectrum  shown  in  Figure  2.17a,  but 
quantitatively  the  effect  of  this  LEBU  is  shown  to  be  much  more  significant.  Immediately 
downstream  of  the  LEBU  trailing  edge,  the  LEBU  modified  deflection  angle  spectrum  is  reduced 
by  approximately  60%  compared  to  the  baseline  deflection  angle  spectrum  below  Sts  =  2. This 
strong  low-frequency  reduction  is  sustained  up  to  approximately  68,  where  the  low-frequency 
(Sts<  2)  end  of  the  deflection  angle  spectrum  begins  to  recover  towards  the  baseline  spectrum. 
Also  just  downstream  of  the  LEBU  device,  there  is  a  slight  20%  increase  in  a  portion  of  the  high- 
frequency  (2  -  78)  end  of  the  deflection  angle  spectrum.  This  increase  disappears  entirely  beyond 
38  downstream  of  the  LEBU  device.  Beyond  3  -  48,  the  entire  high-frequency  range  of  the 
deflection  angle  spectrum  shows  a  reduction  over  the  baseline  of  approximately  40  %  above  Sts  = 
3,  indicating  a  decrease  in  dissipation  relative  to  the  baseline  flow.  This  effect  is  likely  the  result 
of  the  suppression  of  low-frequency  (Sts  <  2)  structures  beginning  immediately  downstream  of  the 
LEBU  trailing  edge  which  yields  an  overall  reduction  in  the  amount  of  energy  cascading  into 
small  scales.  At  x  =  10b,  the  peak  of  the  LEBU  modified  spectrum  is  approaching  Sts  =  1,  which 
is  another  indication  (along  with  the  recovery  of  overall  levels  of  OPDrms  towards  the  baseline 
value)  that  the  aero-optic  reduction  effects  of  this  LEBU  configuration  are  starting  to  wear  off  by 
10& 

Tandem- LEBU  Devices 

Malley  probe  wavefront  measurements  were  also  performed  for  the  tandem  LEBU  devices 
and  overall  levels  of  OPDrms  for  each  tandem  LEBU  device  were  computed  for  streamwise 
locations  x  =  0.8b  to  1 08,  and  the  results,  normalized  by  the  baseline  OPDrms,  are  shown  in  Figure 
2.18.  For  both  tandem  LEBU  configurations,  the  level  of  aero-optic  distortions  just  downstream 
of  the  LEBU  is  reduced  by  approximately  30%.  Further  downstream,  levels  of  OPDrms  for  the  s  = 
48  tandem  LEBU  device  remain  relatively  constant  around  70%  of  the  baseline  value  until  x  ~ 
10b.  Downstream  of  it,  the  aero-optical  distortions  appear  to  be  increasing  slightly  towards  the 
baseline  value.  For  the  s  =  88  tandem  LEBU  device,  OPDrms  continues  to  decrease  further 
downstream  until  reaching  a  region  of  maximum  reduction  between  4-6 8,  where  the  levels  of 
LEBU-modified  OPDrms  are  reduced  by  about  42%  compared  to  the  baseline.  Downstream  of  the 
streamwise  location  of  65,  the  levels  of  OPDrms  begin  to  recover  slightly  towards  the  baseline 
levels,  but  the  rate  of  recovery  is  more  gradual  than  what  was  observed  for  the  /?i  =  0.6b,  hi  =  0.3b 
multi  LEBU  device  (which  had  the  largest  maximum  reduction  in  OPDrms ),  with  the  OPDrms  at 
9- 10b  being  approximately  35%  lower  than  the  baseline  measurement.  This  result  indicates  that 
the  tandem  LEBU  device  results  in  the  most  significant  reduction  in  OPDrms  over  a  longer 
streamwise  extent  than  any  of  the  other  tested  LEBU  configurations. 
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Figure  2.18.  Streamwise  development  of  OPDrms  for  Tandem  LEBU  devices  with  separations  s  = 
4 b  and  8b,  compared  to  the  best-performing  single  and  multi-LEBU  results. 
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Figure  2.19.  Baseline -normalized  spectrum  surfaces  C$  ( x/8 ,  Sts)  for  h  =  0.6  S,  s  =  8b  tandem 
LEBU  devices. 

The  streamwise  evolution  of  the  baseline-normalized  deflection  angle  spectra  are 
presented  in  Figure  2.19  for  h  =  0.6  b,  s  =  8b  tandem  LEBU  device.  The  region  of  suppression  in 
the  low-frequency  end  of  the  spectrum  (Stg<  1)  shows  an  overall  stronger  suppression  of  aero- 
optic  spectrum  compared  to  the  4b  tandem  LEBU  at  the  same  frequency  range  (not  shown). 
Additionally,  for  streamwise  locations  x>  4b  there  is  an  increasingly  strong  reduction  in  the  high- 
frequency  end  of  the  deflection  angle  spectrum.  This  is  likely  is  a  result  a  reduction  in  the  amount 
of  energy  cascading  into  smaller-scales  at  which  turbulence  dissipation  occurs.  For  this  tandem 
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LEBU,  the  baseline-normalized  spectrum  peak  is  also  observed  to  shift  to  lower-frequencies  as 
the  streamwise  location  increases,  moving  from  a  value  around  St$  =  6  to  about  Sts  =  l.At  1  03, 
however,  the  suppression  of  large-scale  structures  ( Sts<  1)  is  still  rather  strong,  so  it  is  possible 
that  levels  of  optical  aberrations  will  remain  low  for  several  more  TBL  thicknesses  before 
beginning  to  relax  to  the  baseline  values. 


Velocity  Measurements 

Earlier  in  this  report  (Section  2.3.1),  velocity  measurements  at  several  streamwise 
locations  in  the  LEBU-modified  TBL  were  presented  along  with  aero-optic  wavefront 
measurements  for  limited  number  of  Single  LEBU  devices.  The  effects  of  LEBUs  on  velocity  and 
wavefront  statistics  were  then  compared,  and  the  results  from  both  measurement  techniques 
indicated  that  frequencies  associated  with  large-scale  turbulent  structures  in  the  TBL  were 
suppressed  by  the  devices.  Both  measurements  also  able  to  identify  LEBU-induced  changes  to  the 
TBL  turbulence  structure  in  the  modified  boundary  layer.  Additional  comparisons  of  velocity  and 
wavefront  statistics  in  the  LEBU  modified  TBL  are  presented  in  this  section  for  selected  LEBU 
configurations  at  a  single  streamwise  location:  x  =  63.  This  streamwise  location  was  chosen  as  it 
is  roughly  in  the  middle  of  the  optical  measurement  section,  and  it  has  been  shown  from  the  data 
presented  in  previous  sections  that  substantial  reductions  in  OPDrms  were  achieved  at  this  location 
for  several  different  LEBU  devices.  The  particular  manipulators  that  will  be  examined  in  this 
section  are  as  follows: 

•  Long  single  LEBU:  /  =  48,  h  =  0.63, 

•  Multi  LEBU:  /  =  1 .63,  ht  =  0 .63,  h2  =  0.3 3,  and 

•  Tandem  LEBU:  /  =  1 .65,  h  =  0.65,  s  =  83. 

Recall  that  the  model  for  calculating  OPDrms  for  compressible  TBLs  have  been  derived 
using  Sutton’s  linking  equation  (1.2)  and  the  Strong  Reynolds  Analogy,  Gordeyev  et  al,  2014. 
The  terms  in  the  SRA  are  easily  arranged  to  show  that  at  a  particular  v-value  prms  in  the  TBL  is 
proportional  to  the  product  U  (y)  urms  (y )  : 

PnJy)  _  r(r  1)iif2  ujyKM .  (2.9) 

Pm  "  Ul 


One  possible  modification  of  the  form  of  the  SRA  given  in  equation  (2.9)  can  be  arrived  at  by 
substituting  the  energy  spectral  density  for  the  root-mean-square  terms  on  both  sides  of  the 
equation.  Since  the  energy  spectral  density  of  any  arbitrary  variable  £  (t)  describes  how  its  total 

variance,  %2ms ,  is  distributed  over  the  component  frequencies  that  c  may  be  decomposed  into.  In 
other  words,  ^ms (/)  =  (/'),  where  S:{f  )  =  F{^2(t)} is  the  one-sided  energy  spectral  density  of 

£  and  F  denotes  the  Fourier  transform.  Therefore,  equation  (2.9)  can  be  re-written  in  the  form 


pL(yJ)=sSRA(yjy 


ul 


where 


SSRA{yJ)  =  F{u{y)u{y,t)]2). 


(2.10) 


(2.11) 


This  decomposes  the  first  term  of  the  integrand  in  the  linking  equation  (1.2)  into  a  function  of 
frequency  and  wall-normal  location.  Assuming  that  the  SRA  is  valid  in  the  LEBU  modified  TBL, 
this  approach  can  be  used  to  identify  where  in  both  frequency  and  wall-normal  space  the  LEBU 
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devices  are  causing  reductions  in  pnns.  To  easily  visualize  this,  the  ratio  of  the  LEBU-modified  to 
baseline  energy  spectra  was  computed: 


C 


SRA 


fy  A 
Ks  s 


i  LEBU 
) SRA 


(y/s, st,) 


(2.12) 


1  BASELINE 
) SRA 


(y/s,stsy 


Here  it  is  emphasized  that  Csra<  1  indicates  a  reduction  in  density  fluctuations  in  the  LEBU- 
modified  boundary  layer,  while  Csra>  1  implies  an  increase  in  density  fluctuations.  Filled  contour 
plots  of  equation  (2. 12)  are  shown  in  Figure  2.20  for  selected  LEBU  devices  at  x/3  =  6. 


st, 

Figure  2.20.  Contour  plot  of  Csra  (y/S,Sts )  at  x  =  63  for  a)  l  =  43  ,h  =  0.63  Long  Single  LEBU,  b) 
hi  =  0.6 3,  hi  =  0.3d  Multi-LEBU  and  c)  .v  =  8d  tandem  LEBU  devices.  LEBU  element  height  is 
marked  by  white  lines,  and  the  solid  black  contour  indicates  where  Csra  (y/3,Sts)  =  1.0. 


The  characteristics  of  the  spectra  surfaces  presented  in  Figure  2.20  are  consistent  with 
other  velocity  and  wavefront  spectra  presented  in  this  and  the  previous  chapters.  The  SRA- 
estimated  fluctuating  density  spectra  show  that  the  LEBU  devices  all  produce  reductions  in  low 
frequencies,  Sts  <  1  (i.e.  large  scales)  below  the  LEBH  height.  Note  that  for  Multi-LEBU  device 
in  Figure  2.20b,  there  is  a  ‘pocket’  of  low-frequency  turbulence  suppression  below  both  LEBU 
elements.  These  spectra  also  show  that  LEBUs  cause  an  increase  in  high-frequencies,  Sts  >  1 
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(corresponding  to  small  scale  structures)  at  wall-normal  locations  y  >  h,  which  corresponds  to 
turbulence  production  due  to  the  increased  mean  velocity  gradient  at  this  location.  The  magnitude 
of  these  increases  is  especially  large  for  long  single  LEBU  and  multi-LEBU  devices. 


Figure  2.21.  Wall-normal  density  correlation  length  functions  for  the  canonical  turbulent 
boundary  layer;  Ai  presented  by  Gilbert  (1982),  A2  measured  by  Rose  &  Johnson  (1982),  and  A3 
computed  from  DNS  of  a  Ree  =  3550  TBL  by  Wang  &  Wang  (2012). 


The  SRA-estimated  prms  (y)  profile  for  the  LEBU-modified  TBL  is  easily  obtained  from 
the  energy  spectrum  Ssra  by  integrating  in  the  frequency  domain: 


Pn 


s(y)= 


/■(/  -  \)p,m: 

u 2 


\sSRA{y,f)df 


,1/2 


(2.13) 


Vo  j 

Note  that  it  can  easily  be  shown  that  equation  (2.13)  is  equivalent  to  equation  (2.9);  however, 
computing  /w  from  equation  (2.13)  allows  for  the  possibility  of  spectral  filtering  if  there  is 
significant  noise  in  the  velocity  measurements. 

It  follows  from  equation  (2.13)  that  the  simplified  form  of  Sutton’s  linking  equation  (1.2) 
can  be  re-written  as: 
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(2.14) 


Remember  here  that  Ap  (y)  is  the  wall-normal  density  correlation  length,  which  can  either  be 
estimated  experimentally  (although  with  some  difficulty)  from  hot-wire  velocity  measurements, 
or  computed  directly  from  CFD  results  (e.g.  Wang  &  Wang,  2012).  Several  different  density 
correlation  lengths  from  the  literature  that  have  been  used  for  modeling  TBL  aero-optic 
aberrations  are  presented  in  Figure  2.21.  These  functions  have  all  been  shown  to  give  similar 
results  for  models  of  canonical  subsonic  and  supersonic  TBLs  (Gordeyev,  et  ah,  2012;  Gordeyev, 
et  al.  2014). 

For  LEBU-manipulated  TBLs,  it  was  not  expected  that  the  functional  forms  of  Ap 
presented  in  Figure  2.21  is  an  appropriate  estimate  for  the  actual  wall-normal  density  correlation 
length,  since  LEBU  devices  significantly  modify  streamwise  and  spanwise  correlations  in 
manipulated  boundary  layers.  However,  if  OPDrms  is  estimated  from  equation  (2.14)  using  the  any 
one  of  the  approximations  for  Ap  for  the  canonical  TBL,  the  influence  of  LEBU-modified  wall- 
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normal  correlation  length  can  be  quantified  by  comparing  predicted  levels  of  reductions  to  those 
measured  using  the  Malley  probe. 

Table  2.3.  Comparison  of  directly  measured  and  indirectly  estimated  values  of  OPD ^  in  LEBIJ 
modified  TBLs  at  x  =  63. 


LEBU  Device  Configuration 

qPDlebu  c 

rms  / 

Measured 

(Malley 

Probe) 

baseline 

From  Eqn.  (2.15),  with  Ap  (y)  = 

Ai  A2  A3 

Single  LEBU  j 

\  1/3=  1.6 
|  h/3  =  0.6 

0.68 

0.82  0.87  0.84 

Long  Single  ! 

LEBU  '  1 

so 

II  II 

0.55 

0.85  0.88  0.86 

Multi-LEBU  < 

'  1/3  =1.6 
h\/3  =  0.6 
h2/3  =  0.3 

V. 

0.64 

0.78  0.75  0.78 

Tandem  LEBU  -< 

"  1/3  =1.6 
h/3  =  0.6 
s/3  =  8.0 

0.55 

0.80  0.83  0.82 

The  reduction  in  OPDrms  relative  to  the  value  predicted  from  baseline  velocity 
measurements  is  given  in  equation  (2.15): 


OPD 


LEBU 

rms 


OPD, 


BASELINE 
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00 
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Mf  ap(tV(t/^) 


(2.15) 


I  o 

where  the  superscripts  ‘ LEBU  and  ‘ BASELINE ’  denote  quantities  for  the  LEBU  modified  and 
baseline  TBLs,  respectively.  Table  2.3  presents  the  results  of  that  the  amount  of  reduction  in 
OPDrms  estimated  from  the  linking  equation  is  significantly  under-predicted  compared  to 
experimental  results.  Depending  on  the  correlation  function  and  LEBU  case,  using  the  canonical 
TBL  correlation  function  over-predicts  OPDrms  in  the  modified  TBL  by  up  to  30%  more  than  the 
actual  value.  This  result  serves  as  a  clear,  although  indirect,  indication  that  LEBU-induced  aero- 
optic  mitigation  is  not  simply  achieved  by  reducing  the  overall  magnitude  of  density  fluctuations 
in  the  TBL.  Rather,  the  de-correlating  ‘plate’  effects  which  were  found  to  be  a  critical  mechanism 
for  LEBU  drag  reduction  (Anders,  1990,  and  references  therein)  are  also  an  important  mechanism 
for  affecting  the  aero-optic  characteristics  of  LEBU-modified  TBLs. 


2.4  Conclusions  and  Discussion 

This  part  of  the  report  presented  experimental  results  of  hot-wire  velocity  and  wavefront 
measurements  of  TBLs,  modified  using  both  single  and  multi-element  LEBU  devices. 
Measurements  of  mean  and  fluctuating  velocity  profiles  downstream  of  LEBU  device,  the  LEBU 
was  found  to  cause  localized  reductions  in  mean  velocity  downstream  of  the  LEBU  device,  as 
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well  as  reductions  in  urms  throughout  the  modified  boundary  layer.  Significant  reductions  in  urms 
were  found  in  the  LEBU  modified  boundary  layer  spectra  -  especially  at  low  frequencies  -  over  a 
significant  portion  of  the  boundary  layer  thickness.  Some  evidence  was  found  that  LEBUs 
introduced  a  significant  amount  of  small-scale  turbulent  structures  into  the  LEBU  wake;  however 
these  are  shown  to  dissipate  quickly  in  the  velocity  spectra.  These  findings  are  consistent  with  the 
results  of  velocity  measurements  from  previous  studies  (Corke,  1981;  Lemay,  et  al.  1990). 

It  was  shown  that  LEBU  devices  are  capable  of  reducing  OPDrms  of  wavefront  aberrations 
caused  by  turbulent  boundary  layers  by  30  -  35%  over  streamwise  distances  of  several  3’s.  The 
LEBUs  caused  this  reduction  by  suppressing  large-scale  turbulent  motions  in  the  modified  TBL. 
These  results  confirm  the  theory  that  careful  manipulation  of  large-scale  turbulent  structures  in 
aero-optically  active  flows  will  be  effective  for  mitigating  aero-optic  aberrations,  and  also 
demonstrate  that  LEBUs  are  an  effective  flow  control  device  for  achieving  this  goal. 

From  careful  analysis  of  wavefront  deflection  angle  spectra,  it  was  also  shown  that  it  is 
possible  to  gain  some  physical  insight  about  the  changes  to  the  modified  TBL.  For  example, 
inspection  of  wavefront  deflection  angle  spectra  showed  that  the  LEBU  devices  suppressed  large- 
scale  turbulent  structures  in  the  manipulated  boundary  layer.  This  result  is  consistent  with  the 
findings  obtained  of  velocity  measurements.  Also,  changes  in  turbulence  production  and 
dissipation,  as  well  as  spectral  features  associated  with  the  growth  of  the  LEBU  wake  were 
identified  through  analysis  of  wavefront  data.  This  supports  the  idea  that  it  is  possible  to  study 
changes  to  compressible  turbulent  flow  structure  using  non-intrusive  aero-optic  wavefront 
sensors,  which  was  initially  put  forth  by  Sutton  (1969). 


x/6 

Figure  2.22.  Baseline  normalized  OPDrms  downstream  of  selected  “best”  LEBU  device 
configurations. 

After  analyzing  all  cases  from  the  parametric  investigation,  the  optimal  tested  single 
LEBU  configuration  was  identified  as  the  l  =  43,  h  =  0 .63,  Tong’  LEBU  device.  This 
configuration  resulted  in  a  maximum  reduction  in  OPDrms  of  approximately  45%,  and  also  gave 
sustained  reductions  over  40%  for  63  in  the  streamwise  direction,  and  reductions  in  OPD rms  over 
30%  for  93.  The  word  ‘tested’  is  emphasized  above,  as  the  optimal  length  found  in  this  parameter 
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space  resided  on  the  edge,  which  leaves  open  the  possibility  that  LEBUs  of  even  longer  chord 
length  will  give  even  larger  reductions  in  OPDrms,  although  LEBUs  with  increasing  chord  length 
will  result  in  additional  device  drag. 

The  effects  of  Multi-LEBU  and  Tandem-LEBU  devices  on  aero-optic  aberrations  in  the 
TBL  were  also  investigated  using  a  number  of  configurations.  A  summary  of  results  for  the  best¬ 
performing  types  of  LEBU  are  presented  in  Figure  2.22.  The  results  demonstrated  that  the  best¬ 
performing  Tandem-LEBU  configuration  for  aero-optic  mitigation  is  two-element  tandem  LEBU 
with  l  =  1 .65,  h  =  0.6 5,  and  5  =  8d.  This  configuration  gave  the  maximum  reduction  in  OPDrms  of 
more  than  40%  downstream  of  the  device  to  about  55,  and  reductions  greater  than  about  35%  over 
a  large  streamwise  range  between  Id  and  1  Or)'.  Similar  reductions  in  OPDrms  of  35%  were  also 
shown  for  a  multi-LEBU  device  with  /  =  1 .65,  h\  =  0.65,  /?2  =  0.3d,  although  the  streamwise 
extent  of  the  reduction  was  limited  to  about  4d. 

Additional  velocity  measurements  were  also  obtained  at  a  single  streamwise 
station  for  several  selected  single,  tandem,  and  multi-LEBU  device  configurations.  The 
fluctuating  density  spectra  were  estimated  for  each  of  these  devices  using  the  Strong  Reynolds 
Analogy,  and  the  results  were  showed  that  significant  reductions  in  the  contributions  from  large 
scales  below  the  LEBU  devices,  and  increased  small  scale  turbulent  fluctuations  above  the  LEBU 
device  that  are  consistent  with  increased  production  on  the  upper  side  of  the  LEBU  wake. 
Estimates  of  OPDrms  were  also  computed  from  velocity  data  through  Sutton’s  linking  equation, 
and  the  results  gave  indirect  evidence  that  changes  in  wall-normal  correlation  length  are  also  a 
key  mechanism  responsible  for  reducing  aero-optic  aberrations  in  the  LEBU-modified  TBL. 

Levels  of  aero-optic  mitigation  for  two  tested  multi-element  LEBU  configurations  were 
shown  to  be  better  than  levels  of  mitigation  achieved  by  the  /  =  1 .65  single-element  LEBU  device. 
However,  the  long  l  =  45  LEBU  performed  better  than  both  of  the  tandem  and  multi  LEBU 
devices,  although  only  just  slightly.  These  results  demonstrate  that  aero-optic  aberrations  caused 
by  the  TBL  can  be  easily  and  significantly  suppressed  by  30  -  40%  over  several  5  in  the 
streamwise  direction  using  Large  Eddy  Break-Up  devices.  The  mechanical  simplicity  of  LEBU 
devices,  along  with  the  experimental  demonstration  of  their  aero-optic  mitigation  abilities  in  this 
work,  makes  them  very  good  candidates  for  incorporation  into  the  design  of  airborne  directed 
energy  systems  that  must  contend  with  aberrations  caused  by  compressible  TBLs. 
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3.  Wall  Cooling  Effects 


3.1  Original  Aero-Optic  Model  for  Wall  Temperature  Effects 


While  aero-optical  effects  caused  by  compressible  boundary  layers  have  been  studied 
since  1950s,  most  of  the  studies  have  not  previously  accounted  for  the  potential  temperature 
mismatch  between  the  underlying  wall  and  freestream  temperature.  So  there  is  a  need  to 
investigate  the  impact  that  heat  transfer  at  the  underlying  wall  of  a  subsonic,  compressible, 
turbulent  boundary  layer  has  on  optical  aberrations. 

Wyckham  and  Smits  (2009)  used  a  high-speed  two-dimensional  Shack-Hartmann  sensor 
to  study  aero-optical  properties  of  subsonic  and  supersonic  boundary  layers.  Although  they  did 
not  perform  heated  wall  experiments,  they  used  a  bulk-flow  analysis  and  developed  a  scaling  law 
for  aero-optical  aberrations  for  non-adiabatic  walls  as  a  function  of  the  freestream  Mach  number, 
the  local  skin  friction  coefficient,  C/,  and  the  ratio  of  the  freestream  static  temperature,  Too,  to  the 


wall  temperature,  Tw,  OPDnns  ~  KGDSMlJcJ 


% 

.  From  this  scaling  law  it  follows  that  aero- 


optical  distortions  should  decrease  when  the  wall  is  heated.  These  predictions  contradicted  the 
preliminary  experimental  studies  conducted  in  Cress  et  al  (2012),  where  it  was  shown  that  aero- 
optical  distortions  increase  when  the  wall  is  heated.  White  and  Visbal  (2012)  have  performed 
Large  Eddy  Simulations  of  aero-optic  aberrations  caused  by  compressible  turbulent  boundary 
layers  for  heated  and  cooled  walls  and  also  have  shown  that  aero-optical  aberrations  increased 
when  the  wall  was  heated.  So,  there  is  a  need  for  a  statistical  model  for  the  OPDnils  as  a  function 
of  subsonic  Mach  number  and  moderate  temperature  difference  between  the  wall  and  the 
freestream. 


In  the  previous  AFORS-funded  work,  the  Extended  Strong  Reynolds  Analogy  was  used  to 
derive  the  model  for  aero-optical  distortions  with  non-adiabatic  wall  was  derived, 


f  f 

=  A(pJPsL)S"Ml  1  +  C, 

>  V 


Cress  (2010)  reported  that  using  velocity  measurements  for  a  boundary  layer  with  M=  0.5,  C\  and 
C2  were  theoretically  predicted  to  be  6.38  and  10.28,  respectively.  A  linearized  form  of  equation 
(3.1)  was  also  found  to  work  well  for  modeling  wavefront  measurements  for  wall  heating,  since 
for  A7>  0,  the  third  term  of  (3.1)  is  much  smaller  compared  than  the  first  two; 


OPDrm,=A(pJpSLyMl  1 +  D 


AT>0. 


where  the  empirical  constant  D\  =  C\I2  (Cress,  2010). 

OPDnns  computed  from  Malley  probe  measurements  of  heated  wall  TBLs  also  showed 
good  agreement  between  experimental  results  and  the  model  equations  (3.1)  and  (3.2),  as  shown 
in  Figure  3.1  for  a  wide  range  of  wall  temperatures  and  Mach  numbers.  These  experiments 
demonstrated  that  OPDrms  could  be  successfully  modeled  using  the  linear  scaling  relationship 
over  a  wide  range  of  wall  temperature  differences  and  subsonic  Mach  numbers.  The  results  of 
these  tests  also  highlighted  the  fact  that  even  moderate  values  of  wall  heating  could  significantly 
increase  wavefront  aberrations  above  levels  encountered  for  adiabatic  wall  temperature  conditions 
only  (Cress,  2010). 
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Figure  3.1.  OPD ^  versus  the  linearized  heated  wall  temperature  scaling,  equation  (3.2),  using 
empirically  determined  D\  constants,  [image  source  Cress  (2010)] 


3.2.  Updated  model  for  non-adiabatic  boundary-layer  aero-optical  distortions 

For  nearly  as  long  as  turbulent  boundary  layers  have  been  studied,  those  with  a  moderate 
heat  transfer  at  the  wall  have  been  the  subject  of  experimental  and  computational  research, 
especially  applied  to  supersonic  and  hypersonic  boundary  layers;  see  Smits  &  Dussauge  (1996), 
Spina  et  al  (1994)  for  summaries  of  the  effects  of  heat  transfer  in  a  compressible,  turbulent 
boundary  layer.  A  significant  advancement  in  the  field  came  in  1962  from  Morkovin  (1962) 
presentation  of  the  Strong  Reynolds  Analogy  (SRA),  which  presumes  that  p ’  is  negligible.  From 
this  analogy  between  the  Reynolds-averaged  form  of  the  energy  and  momentum  equations,  a 
relationship  between  the  fluctuating  static  temperature  and  fluctuating  velocity  can  be  shown  to 
be, 

=  -A(y)[y-\)M 2^~  (3.3) 

where  A(y )  takes  into  account  the  stress  integral  distribution  in  the  boundary  layer  (Smits  & 
Dussauge,  1996);  in  the  original  SRA  A(y )  =  1.  Both  experiments  (Smits  &  Dussauge,  1996; 
Debieve,  1983;  Smith  &  Smits,  1993)  and  DNS  simulations  (Guarini  et  al,  2000;  Duan  et  al, 
2010)  have  shown  that  A(y)  ~  1  for  y/8  <  0.6  and  increases  in  the  outer  layer. 

However,  for  non-adiabatic  wall  conditions,  experimental  (Debieve  et  al  1997)  and 
computational  (Duan  et  al,  2010;  Liu  &  Pletcher,  2007)  studies  have  shown  that  SRA  fails  to 
predict  the  correct  temperature  fluctuations.  Walz  (1969)  proposed  a  form  of  the  enthalpy 

equation  such  that  h  =  h(ll ) ,  where  fluctuations  in  the  total  temperature  are  not  ignored  (Smits 
&  Dussauge,  1996);  and  the  following  relationship  for  mean  static  temperature  was  found, 
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where  T ,  ,  Tw  ?  and  Tr  are  the  Favre-averaged  static,  freestream,  wall,  adiabatic  and  recovery 


temperatures,  respectively,  U  is  the  Favre-averaged  mean  velocity,  is  the  freestream  Mach 
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number  and  r  -  (Tr  -  Tf ) /(T0  —Tx)  is  the  recovery  factor.  In  turbulent  boundary  layers  with  air 

as  the  fluid  in  motion,  r  is  typically  about  0.89.  Equation  (3.4)  is  known  by  several  names:  the 
modified  Crocco  relation,  the  Walz  equation  or  the  Extended  SRA  (ESRA). 

Linearizing  this  equation  in  the  case  of  small  fluctuations  (Debieve  et  al,  1997),  a  relation 
between  temperature  and  velocity  fluctuations  becomes, 


r"  T  -T  ( u"  ^ 

rms  _  r  w  1/1  rms 

T  T  II 

1oo  1  oo  \U  ao  J 


rir-ml -jP 


(3.5) 


where  T”ms  and  u"ms  are  the  Favre-averaged  fluctuating  temperature  and  velocity,  respectively. 

Cebesi  &  Smith  (1974)  derived  a  very  similar  expression  with  A(y )  =  1.  In  general,  equation  (3.5) 
can  be  written  in  the  following  form, 


T  —  T  I  u  |  ,  Uu 

-J— ^  jr  -A(y)(r- ml—f 

*  on  V  ^  nr,  ^  nr. 


Strictly  speaking,  equation  (3.6)  simply  postulates  the  relationship  between  Tvm<,  and  «rms,  while 
ESRA,  in  addition  to  equation  (3.6),  also  requires  that  the  temperature  and  velocity  fluctuations 
are  perfectly  anti-correlated.  This  was  shown  not  to  be  true,  see  for  example  Debieve  et  al  (1997), 
Guarini  et  al  (2000),  Duan  et  al  (2010),  Gaviglio  (1987). 

The  original  model,  equation  (3.2),  was  derived  for  subsonic  boundary  layer  only  and 
assuming  a  particular  relations  between  the  temperature  and  the  velocity,  equation  (3.6).  During 
this  work,  the  model  was  generalized  and  extended  to  a  supersonic  regime.  The  derivation  is 
provided  below. 

From  the  linking  equation  (1.2)  it  follows  that  if  the  density  fluctuations  and  their 
correlation  lengths  across  the  boundary  layer  are  known,  optical  distortions  can  be  calculated. 
From  the  ideal  gas  law,  p  =  pRT,  the  density  fluctuations,  p\  are  related  to  pressure  fluctuations, 
p\  and  temperature  fluctuations,  T .  In  the  case  of  small  fluctuations,  it  can  be  written  as 
p' /  p(y )  =  p' /  p(y) -T' I T (y) .  While  the  pressure  fluctuations  in  adiabatic  boundary  layers  have 


been  shown  to  be  several  times  smaller  than  the  temperature  fluctuations  Wang  &  Wang  (2012), 
Smits  &  Dussauge  (1996),  Spina  et  al  1994),  for  cooled  walls  the  temperature  fluctuations  are 
smaller  than  in  the  adiabatic  case,  Duan  et  al  (2010).  So,  the  pressure  fluctuations  can  be 
comparable  with  the  temperature  fluctuations. 

Using  equations  (3.4)  and  (3.6),  letting  A T=  Tw  -  Tr,  and  replacing  fluctuating  values  with 
root-mean-square  values,  the  following  expression  for  can  be  found, 


2A(y){r-l)Ml 


AT  |  U 


A(y){y-\)M: 


(3.7) 


where  U(y )  is  the  mean  local  streamwise  velocity.  Thus,  the  equation  of  state  can  be  used  to 
compute  the  density  fluctuations,  assuming  no  correlation  between  the  pressure  and  temperature 
(we  will  discuss  this  assumption  later), 


(3.8) 


Here  the  local  pressure  is  assumed  to  be  constant  across  the  boundary  layer.  The  local  density  can 
be  found  from  the  equation  of  state,  p(y)  =  Px  /(RT(y))  =  p^iT^  / T(y)) ,  where  the  mean 


temperature  profile  can  be  computed  from  Equation  (5a), 
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For  non-adiabatic  boundary  layers,  the  fluctuating  velocity  component,  compensated  for  density 


changes  near  the  wall  using  the  Van  Driest  transformation, 


Kms  p(y) 
Ur  V  Pv, 


was  found  to  be  mostly 


unchanged  over  a  wide  range  of  Mach  numbers  (Spina  et  al,  1994,  Guarini  et  al,  2000,  Duan  et  al, 
2010),  so  a  fluctuating  velocity  profile  from  the  adiabatic  boundary  layer,  Guarini  et  al  (2000), 
can  be  used, 

Urmsiy)  \piy)  _  urms(y)  p(y )  _  urms(y)  |  Tm  _  ^  ^  (3  1Q) 
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We  further  assume  that  the  mean  velocity  profile  is  independent  of  the  wall  temperature  as, 

U(y)/Ux=f(y/S )  (3.11) 

Finally,  the  fluctuating  pressure  profile  from  Guarini  et  al  (2000)  will  be  used, 

Prms  /(/VV)  =  prms  /( pMiCf  12 \px  !pw))  =  prms  12))  =  h{y / 8)  (3.12) 

where  the  functions/  g  and  h  are  presented  in  Gordeyev  et  al.  (2015).  Substituting  equation  (3.7) 
into  equation  (3.8)  and  using  equations  (3.10),  (3.11)  and  (3.12)  gives  the  following  relationship 
for  ycws  in  terms  of  the  velocity  and  the  temperature  profiles  in  the  wall  normal  direction  for  a 
given  AT  yields, 
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Finally,  substituting  equation  (3.13)  into  equation  (1.2)  results  in  the  following  relationship, 


OPDrms  =  A0KGDpJJCf  M:  +  C,  —M2  +  C2  — 


(3.14) 


where, 
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Note  that  for  the  adiabatic  wall  boundary  layer,  i.e.,  AT  =  0,  equation  (3.14)  reduces  to  the 
experimentally-proven  scaling  relation,  OPDrms  ~  P^sJcJmI  (Gordeyev  et  al  2012,  2014). 

From  equation  (3.14),  OPD/sis  a  quadratic  function  of  AT  and  reaches  a  minimum  at 
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For  positive  AT)  equation  (3.14)  can  be  rearranged  in  the  following  manner, 


OPD„  =  A0KCDpjJc;  Mi+D,^ 
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where  Dj  =  C\/2  and  D2  =  C3  -  (Ci/2)  .  For  positive  temperature  differences,  the  second  term  in 
the  square  brackets  in  equation  (3.17)  will  be  shown  to  be  much  smaller  than  unity  and  equation 
(3.17)  can  be  further  simplified  as, 


OPD„,„  =  A„KaDp,SyCt 


f  A 

Ml+D,— 

V  j 


Ar>o 


(3.18) 


Several  important  conclusions  can  be  drawn  from  analyzing  the  model  predictions: 


1 .  Cooling  the  wall  to  Armin  should  significantly  reduce  the  aero-optical  distortions. 

2.  Amount  of  cooling,  A rm;n,  is  proportional  to  a  square  of  the  freestream  Mach  number. 

3.  A Tmin  l{T.rM2r  )  primarily  depends  on  the  choice  of  A(y). 

4.  Finally,  heating  the  wall  will  increase  the  aero-optical  distortions. 


In  summary,  the  model  provides  testable  predictions  about  choices  for  A(y)  and  A (y),  and 
comparing  experimental  results  to  the  model  predictions  should  provide  the  evidence  to  support 
(albeit  indirectly)  a  particular  choice  of  both  A(y)  and  A. 

3.2.1  Cooled  Wall  Results 

In  the  previous  AFOSR-funded  investigation,  effects  of  the  wall  heating  were  primarily 
investigated,  with  only  few  experimental  points  for  a  cooled  wall.  In  this  work  the  cooling  effects 
were  thoroughly  documented  and  presented  below. 
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Figure  3.2.  OPDnils  normalized  by  the  OPDmis  value  with  A T=  0  versus  AT the  best-fit 

using  equation  (3.19)  and  the  model  predictions,  equation  (3.14),  using  Ai(y)  and  A( y)  with  and 
without  the  pressure  term. 


For  negative  temperature  differences,  the  full  scaling  relationship,  equation  (3.14)  should  be 
used.  Factoring  out  M2  from  the  right  hand  side  of  equation  (3.14)  gives  the  normalized  OPDmis 


expression, 


OPDms(Ar) 
OPDrms(Ar  =  0) 


1  +  C, 


AT  +  c 

[KMl)  2 


AT 


(3.19) 


As  discussed  earlier,  for  negative  temperature  differences,  aero-optical  distortions  should  be  at 
their  minimum.  The  normalized  OPDmis  data  for  three  Mach  numbers  are  plotted  versus 
AT  /(Ty  M l )  in  Figure  3.2  and  the  strong  effect  of  the  difference  between  the  wall  and  adiabatic 
wall  temperature  is  apparent  in  this  figure.  All  experimental  results  successfully  collapse  onto  one 
curve  and  the  minimum  of  approximately  0.4  in  the  normalized  OPDmis  value  is  clearly  observed 
around  AT /(TyM2yj)  =  -0.4 .  In  other  words,  the  optical  aberrations  in  the  turbulent  boundary 
layer  were  decreased  by  60%  at  this  temperature  difference.  This  is  a  dramatic  decrease  in  the 
magnitude  of  optical  aberrations  and  it  provides  a  promising  passive  way  to  significantly  reduce 
aero-optical  distortions  caused  by  turbulent  boundary  layers.  The  results  also  confirm  the  model 
prediction  that  the  modest  wall  cooling  reduces  OPDnils.  Performing  a  least-square  fit,  the 
constants  Ci  and  C2  were  found  to  be  4.1  and  5.13,  respectively,  equation  (3.19)  with  these 
constants  is  plotted  in  Figure  3.2,  showing  a  fairly  good  prediction  of  the  normalized  OPDrms- 
values  for  the  range  of  AT /(TxMl)  between  zero  and  -0.7;  below  -0.7  experimental  results  are 
higher  than  equation  (3.19)  predicts,  indicating  a  possible  departure  from  some  of  the  assumptions 
used  in  deriving  the  model.  The  value  of  Q  agrees  within  experimental  error  with  the  one 
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obtained  from  the  heated- wall  experiment.  Also,  the  second  term  in  equation  (3.17)  is  always  less 
than  D2  /(2 D?)  =  0.03,  verifying  the  assumption  made  in  deriving  equation  (3.18). 

The  model  requires  the  knowledge  of  the  velocity  profile,  equations  (3.10)  and  (3.11),  the 
pressure  fluctuation  profile,  equation  (3.12),  the  function  A(y),  defined  in  equation  (3.6),  and  the 
density  correlation  length,  A(y).  The  mean  velocity  profile  was  measured  using  the  hot-wire  at  Mx, 
=  0.4,  while  the  fluctuation  velocity  and  pressure  profiles  are  taken  from  the  DNS  simulations, 
Guarini  et  al  (2000).  Three  different  choices  of  the  density  correlation  lengths,  shown  in  Figure 
2.21,  were  used.  For  the  d -function,  three  choices  are  tested:  A(v)  =  r,  A(y )  =  1  and  A(y),  taken 
from  Smits  &  Dussauge  (1996). 


Table  3.1.  Various  parameters  predicted 


by  the  model. 


A(y) 

A(y) 

Ao 

Ci 

c2 

A7;in 

OPD,„JATmm) 

OPDmJATmm) 

OPD,„JAT  =  0) 
,h  =  0 

T,Ml 

OPDrms  (AT  =  0) 

=  r 

Ai 

0.17 

6.41 

10.99 

-0.29 

0.25 

0.10 

=  r 

A2 

0.15 

6.84 

12.69 

-0.27 

0.27 

0.14 

=  r 

A3 

0.13 

6.59 

11.73 

-0.28 

0.27 

0.12 

=  1 

Ai 

0.19 

5.77 

8.80 

-0.33 

0.23 

0.10 

=  1 

A2 

0.17 

6.17 

10.18 

-0.30 

0.27 

0.14 

=  1 

A3 

0.15 

5.94 

9.41 

-0.32 

0.25 

0.12 

A(y) 

Ai 

0.20 

5.46 

7.92 

-0.34 

0.25 

0.14 

A(y) 

A2 

0.18 

5.91 

9.44 

-0.31 

0.27 

0.17 

A(y) 

A3 

0.15 

5.63 

8.54 

-0.33 

0.27 

0.17 

Table  3.1  presents  values  of  A0,  Q,  C2,  ATmia and 
OPDnm(ATmin )  /  OPDnm(AT  =  0) ,  predicted  by  the  model  for  different  selected  A(y)  and  A(y). 
Using  results  from  Table  3.1,  the  closest  prediction  of  the  experimentally-measured  values  of  Ci, 
C2,  the  location  and  the  value  of  the  minimum  OPDnils  is  given  with  Ai  or  A3  and  A( y). 
Theoretical  predictions  with  and  without  the  pressure  term  are  also  shown  in  Figure  3.2.  The 
theoretical  model  predicts  the  location  of  the  minimum  OPDnils  for  A7min  l(TyM2fj )  =  -0.34,  which 
is  within  15%  of  the  experimentally-observed  value.  The  theoretical  model  underpredicts  the 
minimum  value  of  OPDnils  by  30%.  One  possible  source  for  this  discrepancy  is  that  in  deriving 
the  model,  the  cross-term  (T'  p''}  was  not  included.  However,  this  term  is  not  negligible,  as  the 

temperature  fluctuations  are  partially  due  to  the  total  temperature  fluctuations  (Debieve  et  al 
1997;  Guarini  et  al  2000),  which  are,  in  turn,  the  consequence  of  the  pressure  fluctuations.  So  the 
inclusion  of  the  cross-term  would  most  probably  increase  the  level  of  density  fluctuations  and 
would  result  in  higher  values  of  OPDrms. 

The  model  supports  the  choice  of  A(y),  as  all  other  choices  lead  to  even  smaller  values  of 
ATmin /{TXM l) ,  see  Table  3.1.  Notice  that  the  choice  of  A  does  not  significantly  affect 
AT  .  /(t  M2). 

Figure  3.2  also  shows  the  importance  of  including  the  pressure  term  into  the  model,  as 
without  it  the  model  predicts  a  much  lower  value  of  0.17  for  OPDrms(ATmia)  /  OPDrms(AT  =  0) . 


39 


Notice  that  the  inclusion  of  the  pressure  term  is  important  only  in  the  temperature  range  of 
the  lowest  OPDrms  and  does  not  significantly  alter  OPDrmS  for  the  canonical  adiabatic-wall 
boundary  layer.  As  discussed  earlier,  it  is  generally  accepted  that  in  adiabatic-wall  boundary 
layers  the  temperature  fluctuations  are  several  times  larger  than  the  pressure  fluctuations,  so  the 
pressure  fluctuations  can  be  neglected.  But  from  the  presented  model,  equation  (3.7),  it  follows 
that  in  the  moderately-cooled  boundary  layers  the  temperature  fluctuations  are  suppressed,  while 
the  pressure  fluctuations,  which  are  related  to  velocity  fluctuations,  are  not  significantly  modified. 
Thus,  the  cooled-wall  boundary  layer  might  be  a  better  suitable  flow  to  study  the  pressure 
fluctuations  inside  the  boundary  layer. 

A  final  comment  about  the  pressure  term  in  the  model  is  that  its  relative  importance  is 
proportional  to  Cf,  see  equation  (3.16),  and  it  would  be  even  more  important  at  low  Reynolds 
numbers. 

As  mentioned  before,  the  only  other  model  for  including  temperature  effects  to  predict 
aero-optical  distortions  of  a  boundary-layer  is  the  model  proposed  in  Wyckham  &  Smits  (2009). 
From  their  model  it  follows  that  optical  aberrations  are  inversely  related  to  the  wall  temperature; 
their  model  predicts  that  the  value  of  OPDrms  will  decrease  as  the  wall  temperature  is  increased. 
Clearly,  the  presented  experimental  data  pointedly  contradicts  this  result.  A  close  inspection 
shows  that  the  model  from  Wyckham  &  Smits  (2009)  is  based  upon  the  SRA,  not  the  “extended” 
SRA,  and  assumes  that  total  enthalpy  is  constant  throughout  the  boundary  layer,  and  therefore 
does  not  allow  the  total  temperature  to  vary. 


3.2.2  Convective  speeds 


The  presented  model,  equation  (3.13),  estimates  the  boundary-layer  density  profile  for 
different  wall  temperatures.  One  can  make  a  plausible  assumption  that  the  density  structure 
convects  at  the  local  velocity  and  the  observed  overall  convective  speed  of  the  aero-optical 
structure  is  related  to  the  velocity  integral  weighted  by  the  density  field.  So,  we  can  speculate  that 
the  measured  convective  speed  of  aero-optical  structures  can  be  presented  as, 


uc  _  \\PrmXy^T)\Kyld)dy 

u *  "  \\prms(yAT)\dy 


(3.20) 


It  is  worth  noting  that  the  similar  approach  was  used  to  correctly  predict  the  experimentally- 
observed  increase  in  the  convective  speed  of  aero-optical  structures  in  a  boundary  layer  at 
supersonic  speeds,  Gordeyev  et  al  (2012),  Gordeyev  et  al  (2015).  The  model  prediction,  equation 
(3.20),  with  A( y),  is  plotted  in  Figure  3.3  as  a  function  of  the  wall  temperature,  along  with 
experimental  measurements  of  the  convective  speeds.  While  the  particular  choice  of  the  integral 
in  equation  (3.20)  is  admittedly  somewhat  ad  hoc,  the  model  properly  predicts  all  the 
experimentally-observed  trends,  including  the  reduction  of  the  convective  speed  for  large  negative 
A  T!(t„MI)<  -0.5  and  even  absolute  values  of  the  convective  speed  for  the  adiabatic  wall 
boundary  layer.  From  the  model  one  can  see  that  for  large  negative  temperature  differences,  the 
mean  temperature  profile  near  the  wall,  equation  (3.9),  was  modified,  compared  to  the  adiabatic 
wall  case.  This  led  to  the  increase  of  the  density  fluctuations  near  the  wall,  while  in  the  outer  part 
of  the  boundary  layer  the  density  fluctuations  were  still  suppressed.  Therefore,  the  slower-moving 
structures  near  the  wall  were  optically  “amplified”,  leading  to  the  decrease  in  the  overall 
convective  speed  of  aero-optical  structure.  The  model  prediction  without  the  inclusion  of  the 
pressure  term  is  also  shown  in  Figure  3.3  and  shows  a  lesser  agreement  with  the  data.  Thus,  it  also 


40 


indicates  that  the  pressure  term  has  to  be  included  to  study  density  fluctuations  in  moderately- 
cooled  boundary  layers. 
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Figure  3.3.  Convective  velocity  of  aero-optic  aberrations  for  full  wall  cooling  as  a  function  of 
wall  temperature  and  comparison  with  the  theoretical  predictions. 


3.3  Extension  to  supersonic  speeds 

The  model  reasonably  predicts  all  the  important  aero-optical  characteristics  of  non- 
adiabatic-wall  subsonic  boundary  layers;  it  also  correctly  predicts  OPDrms  for  supersonic  adiabatic 
boundary  layers,  Gordeyev  et  al  (2012).  Because  the  density  correlation  lengths  (Gao  et  al,  2013), 
and  normalized  pressure  variations  (Duan  &  Choudhari,  2013),  are  similar  for  subsonic  and 
supersonic  speeds,  the  model  predictions  can  be  extended  to  the  supersonic  non-adiabatic  regime. 
Results  of  the  predicted  levels  of  OPDrms  at  different  Mach  numbers  and  temperature  differences, 
adjusted  to  match  experimental  results  at  subsonic  speeds  are  shown  in  Figure  3.4.  At  supersonic 
speeds,  the  OPDrms-level  becomes  a  non-linear  function  of  the  wall  temperature  for  the  heated 
boundary  layer  and  the  relative  increase  in  OPDrms  becomes  less  pronounced  for  Mach  numbers 
larger  than  one.  OPDrms  still  has  a  minimum  value  around  AT  /(tx  MI  )=  -0.4  for  all  Mach 
numbers;  also  the  relative  reduction  in  OPDrms  is  similar  for  both  subsonic  and  supersonic 
boundary  layers.  Note  that  for  M^  =  2.0  and  M^  =  3.0  the  lines  in  Figure  10  do  not  extend  to  large 
negative  temperatures  because  the  wall  temperature  eventually  reaches  the  absolute  zero  and 
cannot  be  cooled  further. 

Overall,  the  wall  cooling  still  should  be  an  effective  way  to  reduce  aero-optical  distortions 
of  the  boundary  layer  at  low  supersonic  speeds,  M*,  <  3.  The  optimal  cooling  is  expected  to  be 
around  AT  /(TxMl)=  -0.4  with  the  reduction  in  OPDrmS  by  a  factor  of  two  or  so. 


3.4  Results  for  Partial  Wall  Cooling 

The  predictions  from  the  model  raise  the  possibility  of  using  wall  cooling  as  a  method  for 
aero-optic  mitigation.  In  practical  applications  where  this  is  of  potential  use,  however,  the  energy 
expense  of  cooling  of  the  full  TBL  development  length  in  an  airborne  application  might  be 
impractical  (Cress,  2010).  One  alternate  means  of  applying  wall  cooling  flow  control  might  be  to 
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apply  cooling  over  only  partial  lengths  of  the  TBL  development  length  upstream  of  the  optical 
aperture. 


Adjusted  Model 


Figure  3.4.  Adjusted  model  predictions  (with  Ai  and  A(y))  at  subsonic  and  supersonic  speeds. 

In  addition  to  the  full  wall  cooling  case,  several  different  types  of  partial  wall  cooling 
configurations  were  evaluated  using  the  Malley  probe  wavefront  sensor.  A  set  of  near-aperture 
partial  wall-cooling  wavefront  measurements  were  obtained,  in  which  the  cooling  end  location  x\ 
was  fixed  just  upstream  of  the  optical  measurement  section  at  x  =  120  cm,  and  the  location  where 
cooling  began,  xo,  was  varied  in  order  to  change  the  length  of  cooling,  Lcooi  =  x\-  xo.  The  second 
category  of  partial  wall-cooling  tests  was  a  set  of  wavefront  measurements  for  far-upstream 
cooling,  in  which  xo  was  fixed,  and  x\  was  varied  in  order  to  change  the  cooling  length.  The  full 
list  of  the  wall-cooling  cases  tested  in  this  experiment  is  given  in  Table  3.2  along  with  the 
streamwise  locations  of  the  start,  xo,  and  end,  xi,  of  the  cooled  wall  segments.  The  ratio  of  Lcooi  to 
the  TBL  development  length  to  the  Malley  probe  location,  Rcooi  =  LcooilxMp,  are  also  given  in 
Table  3.2. 

Table  3.2,  Parameters  for  full  and  partial  wall  cooling  cases. _ 


Sections 

Cooled 

Case  Description 

Moo 

L COOL 

L DEV 

1,2,3 

Full  Wall  Cooling 

o 

o 

96% 

1 

Far-Upstream  (Short) 

0.35..  0.4 

32% 

2 

Upstream  Strip  (Short) 

0.35..  0.4 

32% 

3 

Near-Aperture  (Short) 

0.35..  0.4 

32% 

1,2 

Far-Upstream  (Long) 

0.4 

64% 

2,3 

Near- Aperture  (Long) 

0.35..  0.4 

64% 
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OPDrms  was  computed  from  partial  wall  cooling  wavefront  measurements  in  the  same 
manner  as  for  full  wall  cooling;  the  results  are  presented  in  Figure  3.5  along  with  the  M=  0.4,  full 
wall  cooling  OPDrms  data  for  comparison.  The  results  for  the  three  longer  partial  wall  cooling 
configurations,  where  Rcooi  =  64%,  are  shown  in  Figure  3.5a.  For  the  far-upstream  64%  cooling 

configuration,  the  minimum  value  of  OPDrms  was  about  0.53OPDrms[AT  =  0),  and  value  of 

AT optimai  was  approximately  -§AlTfM  For  the  same  length  of  wall  cooling  near  the  aperture,  the 
optimal  wall  temperature  was  found  to  be  about  -OAOTJM2,  and  effectiveness  for  aero-optic 
mitigation  is  slightly  better  than  for  the  far-upstream  cooling,  at  0A5OPDrms(AT=  0).  The  optimal 
cooling  temperature  for  the  far-upstream  cooling  case  was  about  20%  lower  than  the  value  found 
for  full  wall  cooling;  however  for  near  aperture  cooling  A  T0ptimai  was  not  altered  significantly. 
Recall  from  the  description  in  Table  3.2  that  the  third  partial  wall  cooling  configuration  was  a 
hybrid  of  short  near-aperture  and  far-upstream  cooling,  each  with  Rco„i  =  32%.  Combined,  this 
results  in  a  total  non-contiguous  cooling  area,  Rcooi,  equal  to  64%  of  the  total  TBL  development 
length.  The  near/far  combined  configuration  yielded  an  optimal  cooling  temperature  estimated  to 
be  in  the  neighborhood  of  A T optimal  =  -OAQTyM2  (although  wavefront  data  was  not  obtained  for 
AT/(TooM2)  =  -0.5  to  -0.3).  For  the  whole  range  of  wall  temperatures  measured  for  this  case, 
OPDrms  is  slightly  higher  compared  to  the  full  wall  cooling  case. 

Data  for  shorter  lengths  of  wall  cooling  are  presented  in  Figure  3.5b  for  three  additional 
partial  wall  configurations:  near-aperture  cooling,  far-upstream  cooling,  and  ‘mid-range’  cooling, 
in  which  the  middle  one-third  of  the  wall  in  the  middle  of  the  TBL  development  test  section  was 
cooled.  For  the  short  length  of  far-upstream  cooling,  A T0ptimai  was  found  to  be  -0.537CM2,  and 
OPDrms  was  only  reduced  by  approximately  28%  at  this  temperature.  For  the  short  length  of  near 
aperture  partial  cooling,  A T0ptimai  was  about  -0.387TTT2,  where  aero-optic  aberrations  were  only 
reduced  by  about  23%.  These  data  show  that  the  aero-optic  mitigation  performance  is  slightly 
better  for  the  near-aperture  configuration  than  the  far-upstream  case  of  the  same  length,  but  the 
differences  are  not  very  large.  Cooling  a  strip  of  the  wall  in  the  middle  of  the  TBL  development 
section  (i.e.  mid-range  strip  cooling),  resulted  in  a  shift  in  A T optimal  to  a  lower  temperature  (~ 
-0.57’ooAr)  than  was  found  for  full  wall  cooling.  However,  this  configuration  was  found  to  result 
in  a  slightly  larger  reduction  in  OPDrms  (about  a  44%  decrease)  than  either  of  the  other  short 
length  cooling  configurations. 

While  OPDrms  measured  for  a  variety  of  different  partial  wall  cooling  configurations 
possessed  a  functional  form  similar  to  the  one  that  has  been  both  observed  and  predicted  for  full 
wall  cooling,  changes  in  optimal  temperature  value  and  OPDrms  reduction  were  found  to  be  highly 
dependent  on  the  partial  wall  cooling  configuration.  It  was  shown  before  that  both  the  optimal 
temperature  value  and  the  amount  of  aero-optic  mitigation  are  a  function  of  the  model  constants 
Ci  and  C2.  Therefore,  it  is  likely  that  an  effective  means  of  modifying  the  model  for  partial  wall 
cooling  is  to  change  the  model  constant  definitions  so  that  they  are  functionally  dependent  on  the 
length  and  location  of  wall  cooling. 
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A  TKT^M2) 

Figure  3.5.  Normalized  measurements  of  OPDrms  as  a  function  of  wall  temperature  AT/  ( TXjM2) 
for  partial  wall  cooling  cases:  a)  long  lengths  of  cooling  (64%)  and  b)  shorter  lengths  of  cooling 
(32%). 


Convection  velocity  was  also  computed  from  wavefront  measurements,  and  the  results  for 
the  shortest  strips  (Rcooi  =  32%)  for  both  near  aperture  and  far-upstream  cooling  are  plotted  in 
Figure  3.6.  The  results  of  measured  Uc  from  full-wall  cooling  measurements  (from  Figure  3.3) 
are  also  shown  for  comparison.  For  far-upstream  wall  cooling,  the  convection  velocity  is  reduced 
below  the  value  of  0.83  {%,  measured  for  adiabatic  wall  TBLs,  with  Uc  decreasing  monotonically 
with  wall  temperature  to  about  0.736%  at  AT  =  -0. 8 TJV[2' Variation  in  Uc  with  wall  temperature 
was  also  found  in  full  wall  cooling  wavefront  measurements,  however  for  far-upstream  partial 
wall  cooling  the  measured  values  of  Uc  showed  were  consistently  lower  than  those  found  for  full- 
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wall  cooling.  For  near-aperture  cooling,  convective  velocity  is  actually  increased  over  the  value  of 
0.83  U«,  measured  for  adiabatic  wall  TBLs.  As  wall  temperature  was  decreased  to  around 
-0.45  FooM,  Uc  increased  to  a  maximum  value  of  about  0.881/00,  and  as  temperature  was 
decreased  further,  Uc  was  found  to  decrease,  reaching  a  value  of  about  0.78  f/o  around  AT  = 
-0.75  JooM2' 

In  Figure  3.6,  it  is  clear  that  the  functional  dependence  of  Uc  on  AT  for  partial  wall 
cooling  temperature  is  distinctly  different  than  the  functional  dependence  measured  for  full  wall 
cooling.  Measured  values  of  Uc  for  near  aperture  partial  wall  cooling  were  found  to  be  generally 
greater  than  for  full  wall  cooling,  while  Uc  measured  for  far-upstream  cooling  was  generally  less 
than  the  full  wall  cooling  results.  One  possible  explanation  for  these  differences  was  discussed  in 
Chapter  4,  namely  that  step  changes  in  wall  temperature  create  distinct  thermal  sub-layers  within 
the  turbulent  boundary  layer.  Assuming  that  this  is  the  case,  far-upstream  cooling  could  generate 
a  thermal  sub-layer  that  would  affect  density  fluctuations  in  the  outer  part  of  the  turbulent 
boundary  layer.  Conversely,  a  thermal  sub-layer  generated  by  near  aperture  cooling  would  tend  to 
affect  density  near  to  the  wall,  reducing  the  aero-optic  signature  of  this  region.  The 
experimentally  observed  effect  on  Uc  is  consistent  with  this  hypothesis. 


Figure  3.6.  Convective  velocity  computed  from  partial  wall  cooling  wavefront  measurements. 

3.5  Partial  Wall  Cooling  Model 

Based  on  characteristics  observed  in  partial  wall  heating/cooling  literature  and  in 
convective  velocity  measurements  from  different  partial  wall  cooling  cases,  it  is  reasonable  to 
assume  that  the  variation  in  the  effect  of  the  length  and  location  of  wall  cooling  on  aero-optic 
aberrations  is  caused  by  the  differences  between  thermal  sub-layers  generated  by  each  unique 
configuration.  Since  the  model  constants  C\  and  C2  account  for  the  effect  of  wall-temperature  in 
the  TBL  aero-optic  model,  equation  (3.14),  it  was  assumed  that  these  empirical  constants  would 
retain  the  functional  form  found  in  equation  (3.14)  for  partial  wall  cooling,  but  with  C\  and  C2 
being  dependent  on  the  start  location,  x0,  and  end  location  x\,  of  wall  cooling: 
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This  notion  is  supported  by  the  general  characteristics  of  OPDrms  ( AT)  observed  in  experiments 
for  a  variety  of  different  partial  wall  cooling  configurations,  as  shown  in  Figure  3.5. 


Table  3.3  Empirical  constants  computed  for  full  and  partial  wall  cooling  experiments 


Configuration 

Rcooi 

Ci 

c2 

AT 

1  Optimal 

LMl 

opdJat^) 

opd,Jat  =  o) 

Full  Wall 

96% 

4.42 

5.71 

-0.39 

0.38 

Far-Upstream  (Short) 

32% 

2.09 

2.18 

-0.48 

0.71 

Far-Upstream  (Long) 

64% 

3.07 

3.32 

-0.46 

0.54 

Near  Aperture  (Short) 

32% 

2.43 

3.67 

-0.33 

0.77 

Near  Aperture  (Long) 

64% 

3.95 

4.94 

-0.40 

0.46 

Mid-Range  (Short) 

32% 

2.74 

2.74 

-0.50 

0.56 

Near/Far  Combined 

64%* 

3.97 

5.29 

-0.38 

0.51 

* Length  of  cooling  is  not  continuous. 


To  empirically  determine  the  values  of  C\  and  C2  that  corresponded  to  each  partial  wall 
cooling  case  from  experiments,  equation  (3.21)  was  fit  to  the  results  of  each  experiment  over  the 

range  —  0.7  <  AT/fM1)  <  0,  and  the  results  are  presented  in  Table  3.3  along  with  the  values  of 

AToptimai /(Nm2)  and  OPDrms(AT0ptimal)/OPDrms(AT  =  0)  computed  from  the  empirically 

determined  C\  and  C2.For  increasing  values  of  Rc00i,  the  experimentally  determined  values  of  both 
constants  also  monotonically  increased  towards  the  values  obtained  for  full  wall  cooling 
wavefront  measurements.  The  values  of  C\  and  C2  calculated  for  far-upstream  cooling  were  lower 
than  those  found  for  all  other  configurations.  Out  of  all  the  configurations  where  Rcooi  =  32%,  the 
mid-range  case  (where  Section  2  only  was  cooled)  resulted  in  the  largest  aero-optic  mitigation. 
For  cases  where  Rcooi  =  64%,  the  near-aperture  cooling  case  gave  the  biggest  reduction  in  OPDrms. 
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Aberrated 


Figure  3.7.  Schematic  of  a  simple  thermal  sub-layer  used  for  modeling  the  effects  of  partial  wall 
cooling  on  wavefront  aberrations  in  the  TBL. 


Previous  characterizations  of  thermal  sub-layers  caused  by  step  changes  in  wall 
temperature  (Gran,  et  al.,  1974;  Antonia,  et  al.,  1977)  were  used  to  inform  the  modification  of  the 
existing  wall  cooling  model  for  the  effect  of  partial  wall  cooling.  To  simulate  the  effect  of  thermal 
sub-layer  on  boundary  layer  aero-optic  statistics,  the  scalar  AT  in  equation  (3.7)  was  replaced  with 
an  arbitrary  temperature  difference  function  Fat,  which  is  a  function  of  y  and  the  upper  and  lower 
bounds  of  the  thermal  sub-layer,  do  and  Sy. 
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Comparing  equation  (3.22)  to  (3.7)  w  it  becomes  clear  that  Fat  (do,  di,  y)  can  be  thought  of  as  the 
effective  value  of  AT  at  any  point  y  in  the  TBL  profile  for  a  given  wall  cooling  configuration.  It 
follows  then  that  Fat  can  be  written  as  the  product  of  the  wall  temperature  difference  AT  and  a 
configuration- specific  sub-layer  mask  function,  117’  (do,  Si,  v)  :  Fat  (do,  Si,  v)  =  AT^Ut  (do,  Si,  v). 
Note  that  fir  is  bounded  between  0  and  1  and  has  the  same  functional  form  as  Fat  (see  inset  of 
Figure  3.7).  Substituting  this  relation  into  (3.22),  the  temperature-velocity  relationship  for  partial 
wall  cooling  becomes 
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ATxU^S^y) 
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(3.23) 


In  effect,  the  addition  of  a  mask  function  to  the  model  limits  the  application  of  the  ESRA 
temperature-velocity  relation  to  the  region  of  the  TBL  within  the  thermal  sub-layer;  outside  of  the 
sub-layer,  where  IIr=  0,  equation  (3.23)  reduces  to  the  SRA. 

This  temperature- velocity  relation  was  combined  with  equation  (3.8)  to  find  that 
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where  T  (y)  /7T  is  given  in  equation(8a).  Equation  (3.24)  was  then  integrated  using  Sutton’s 
linking  equation  (1.2)  to  calculated  OPDrms.  The  resulting  aero-optic  model  equation  reduces  to 
the  same  functional  form  found  in  equation  (3.14)  with  the  sub-layer  mask  function  being 
grouped  into  the  integrals  for  C\  and  C2. 

ci(*o>*i)=  2^2  ^ J nr , <?i  )A(y)f(y)g2 A(y)dy  (3-25a) 


(3.25a) 


C2(x0,Xi)—  Jn T(y,S0,Sl)g  A {y)dy 


(3.25b) 


where  So  and  Si  are  functions  ofxo  andxi,  respectively. 

A  piecewise  step  function  was  used  to  define  the  sub-layer  mask  function: 

0,  if  y<Sx, 

nr0',5o.^,)=  1. 


UT(y,S0,Sl)=  jl,  if  S{<y<S0,  (3.26) 

0,  if  y  >  S() . 

The  thermal  sub-layer  limits  So  and  S\  were  calculated  at  the  Malley  probe  measurement  location, 
xmp,  using  the  following  approximations: 


S0  =  S(x, 


S]  =  S(* 


XMP  X0 


XMP  X\ 


(3.27) 


These  growth  rates  for  these  internal  sub-layer  boundaries  were  estimated  from  the  results  of  prior 
investigations  of  thermal  step  changes  in  the  wall  boundary  conditions  presented  by  Gran,  et  al. 
(1974),  Antonia,  et  al.,  (1977),  and  others. 

It  follows  from  equation  (3.21)  that  the  optimal  temperature  difference  1ST optimal  at  which 
maximum  aero-optic  mitigation  is  achieved  for  partial  wall  cooling  is  computed 

A  Poptinud  (*0  Ai )  _  1  C\  (xp ,  %i )  ('3  28') 

T.Ml  “  2C,(x0,xJ 

Subsequently,  the  value  of  the  baseline  normalized  OPD at  the  optimal  temperature  is 


npnN0RM (\T  \-  1  1  Q  (Wi) 
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Since  using  the  function  of  A  (y)  from  Smits  &  Dussauge  (1996)  gave  model  predictions 
that  were  closest  to  experimental  results  for  the  full  wall  cooling  model,  this  term  will  be  used  in 
the  partial  wall  cooling  model  predictions. 


□  Far-Upstream  Cooling 
•  Near  Aperture  Cooling 
—  Far:  Ai 

- Near:  Ai 

— -  Far:  A2 

- Near:  A2 

Far:  A3 

. Near:  A3 

C)  ATol>l!mai/(T,MlX  and  d) 


2 

The  effect  of  Ap  (y)  on  model  predictions  of  C\,  Cj,  AT 'optimal! (TJ\T ),  and 
OPDZRM{AT0ptlJ  is  shown  in  Figure  3.8  along  with,  the  values  calculated  from  experimental 


results.  The  choice  of  wall-normal  density  correlation  function,  A p(y),  was  found  to  have  a  strong 
influence  on  the  shape  of  the  functional  dependence  of  model  predictions  on  Rcooi-  All  of  the 
model  predictions  tended  to  over-predict  the  near  aperture  cooling  constants  with  respect  to  the 
experimentally  obtained  values.  The  results  of  the  model  predictions  using  A2  over  predicted 
experimental  results  by  as  much  as  40%,  while  model  predictions  using  Ai  and  A3  only  gave  a 
23%  over  prediction.  For  the  far-upstream  case,  all  three  choices  for  Ap  gave  decent  model 
predictions  of  C\  and  C2  around  Rcooi  =  0.64,  but  significantly  under  predicts  for  Rcoo/  <  0.64,  and 
over  predicts  for  Rcooi  ~  1.  For  all  three  choices  of  density  correlation  length,  the  partial  wall 
cooling  model  gave  estimates  of  A  T0ptimai  or  OPDrms  that  were  not  in  great  quantitative  agreement 
with  experimental  results,  but  had  functional  forms  that  were  consistent  with  trends  observed  in 
experiments. 
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□  Far-Upstream  Cooling 
•  Near  Aperture  Cooling 
—  Far:  Ai 

- Near:  Ai 

— -  Far:  A2 

- Near:  A2 

Far:  A3 
. Near:  A3 


Figure  3.9.  Partial  wall  cooling  model  predictions  for  different  choices  of  Ap  (y),  scaled  to  match 
empirically  determined  full-wall  constants. 


It  is  apparent  that  the  simple  analytical  model  proposed  in  equation  (3.25)  does  a  good  job 
of  qualitatively  modeling  the  effect  of  partial  wall  cooling  on  aero-optic  characteristics. 
Quantitatively,  however,  the  model  has  been  shown  to  be  sensitive  to  different  choices  for  input 
parameters,  and  to  have  generally  poor  agreement  with  experimental  results.  To  correct  this 
quantitative  discrepancy,  the  model  predicted  constants  C\  (xo,  *i)  and  C2  (xo,  xQ  were  scaled  so 
that  they  are  equal  to  the  experimentally  obtained  values  of  Ci  and  C2  for  the  full  wall  cooling 
case  {Rcool  ~  1).  The  results  of  this  scaling  are  shown  in  Figure  3.9,  for  different  values  of  Ap  (y) 
using  the  same  parameters  for  the  calculations  presented  in  Figure  3.8.  The  results  show  that  the 
model  predictions  for  Ai  and  A3  match  experimental  data  well  for  near-aperture  cooling.  For  the 
case  of  far-upstream  cooling,  Ai  performs  the  best,  but  still  under-predicts  C\  and  C2  and  over¬ 
predicts  OPDrms  {AT Optimal)  with  respect  to  the  experimentally  determined  values.  Despite  this 
discrepancy,  the  model  still  does  a  reasonably  good  job  of  predicting  the  A T0ptimai  for  far- 
upstream  cooling. 

This  good  qualitative  agreement,  and  the  prior  literature  on  step  changes  in  wall 
temperature,  indicates  that  the  thermal  sub-layer  assumption  for  modeling  partial  wall  cooling  is 
the  appropriate  approach,  but  the  quantitative  disagreement  demonstrates  that  particular 
methodology  used  in  the  present  study  does  not  fully  capture  the  physics  of  all  partial  wall 
cooling  configurations. 
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Figure  3.10.  Convection  velocity  of  aero-optic  aberrations,  measured  from  partial  wall  cooling 
experiments  and  predicted  from  the  model  equation  (3.30). 

Earlier  in  this  chapter,  it  was  shown  that  equation  (3.24)  predicts  profiles  of  /w  in  the 
non-adiabatic  TBL.  Similar  to  equation  (3.20),  the  convective  velocity,  Uc,  of  wavefront 
aberrations  in  the  partially-cooled  wall  TBLs  can  be  estimated  using  relation 


)  _  \\Prms(y^T’xo ^)\f{y I S)dy 

\\prms(y^T’x^xx)\dy 


(3.30) 


U, 


oo 


The  model-predicted  dependence  of  Uc  on  AT  was  computed  from  equation  (3.30)  for  the  shorter 
(32%)  lengths  of  wall  cooling  for  which  wavefront  measured  values  of  Uc  were  reported  in 
Figure  3.6,  and  the  predictions  are  compared  to  experimental  data  in  Figure  3.10.  Similarly  to 
other  model  predictions,  the  results  shown  in  this  figure  demonstrate  that  while  the  model- 
predicted  values  of  Uc  are  in  poor  quantitative  agreement  with  experimental  results,  the 
configuration-dependent  trends  predicted  by  the  model  are  very  similar  to  those  found 
experimentally. 

3.6  Conclusions  and  Discussion 

The  primary  objective  of  this  work  is  to  present  a  theoretical  model  for  the  fluctuating  density 
field  in  adiabatic,  cooled  and  heated  boundary  layers  based  on  the  Extended  Strong  Reynold 
Analogy.  The  fluctuating  density  field  in  turn  can  be  used  to  predict  the  aero-optical  distortions 
in  adiabatic-  and  non-adiabatic-wall  boundary  layers.  Strictly  speaking,  the  model  presented  here 
does  not  require  a  perfect  anti-correlation  between  the  velocity  and  the  temperature  fluctuations 
but  rather  simply  assumes  a  functional  dependence  between  them.  The  experimental  data  were 
collected  specifically  to  validate  the  theoretical  model;  these  data  were  used  to  test  predictions  for 
heated-  and  cooled-wall  experiments  at  several  subsonic  Mach  numbers.  The  wall  temperature 
effect  on  OPDrms  was  successfully  predicted  by  the  model.  For  a  heated  wall  the  data  validated  the 
model’s  prediction  that  the  level  of  aero-optical  distortion  is  a  linear  function  of  AT /(r^M2 ).  For 
the  cooled  wall  boundary  layer,  aero-optical  distortions  initially  decrease  with  a  negative 
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AT /(Ty Ml),  reach  a  minimum  value  at  the  optical  wall  temperature  of  AT l{TYM2Y)=  -0.4  and 
finally  increase  for  larger  negative  temperature  differences;  this  behavior  was  also  successfully 
predicted  by  the  model.  By  assuming  different  relations  between  the  velocity  and  the  temperature 
fields  and  comparing  the  theoretical  predictions  with  the  experimental  results,  the  particular 
functional  dependence  between  the  velocity  and  the  temperature  fields  was  generally  supported. 
In  addition,  the  comparison  demonstrated  that  the  pressure  term  should  be  included  to  explain  the 
density  fluctuations  for  the  moderately-cooled  boundary  layers.  The  model  also  correctly  predicts 
the  variation  in  the  aero-optical  convective  speed  for  cooled  walls.  Finally,  the  model  predictions 
for  different  temperature  differences  were  extended  to  supersonic  speeds. 

Aero-optical  distortions  are  directly  linked  to  the  density  structure  by  the  linking  equation, 
equation  (1.2),  so  studying  these  aero-optical  distortions  gives  many  important  insights  about 
properties  of  the  density  field  in  boundary  layers,  which  is  a  hard  quantity  to  measure 
experimentally.  Also,  by  varying  the  wall  temperature,  the  temperature  field  can  be  modified 
independently  from  the  pressure  field  and  their  respective  effects  on  the  density  field  can  be 
accurately  and  non-intrusively  measured. 

The  results  show  that  a  heat  addition  or  subtraction  at  the  wall  in  a  subsonic,  compressible 
turbulent  boundary  layer  has  an  important  impact  on  the  optical  aberration  imprinted  on  the 
wavefront  of  a  laser  propagated  through  the  boundary  layer.  Particularly,  it  implies  that  the 
aircraft  skin  upstream  of  the  optical  aperture  should  be  kept  close  to  or  below  the  recovery 
temperature;  otherwise,  any  airborne  system  which  requires  optical  propagation  through  the 
boundary  layer  will  experience  larger  values  of  aero-optical  distortions.  Using  cooled  walls,  aero- 
optical  distortions  might  be  decreased  by  as  much  as  60%,  providing  potentially  important 
practical  means  to  passively  mitigate  aero-optical  effects  in  boundary  layers. 

From  the  deflection  angle  spectra,  it  is  apparent  that  for  the  range  of  AT  values  tested,  a 
moderate  AT  from  the  heated  wall  simply  amplifies  the  density  fluctuations  in  the  boundary  layer 
with  little  or  no  effect  on  the  velocity  statistics  of  the  boundary  layer.  Further,  the  amplitude  in  the 
deflection  angle  spectrum  varied  linearly  with  the  temperature  difference.  Thus,  increasing  the 
wall  temperature  will  not  alter  the  aero-optical  structure,  but  rather  will  simply  amplify  it 
optically,  at  least  in  the  temperature  range  tested.  These  results  suggest  that  low  speed  flows  with 
very  low  levels  of  optical  aberrations  can  be  made  optically  active  by  introducing  a  temperature 
mismatch  between  the  wall  and  freestream;  therefore,  optical  aberrations  that  are  unobservable  at 
AT  =  0  can  be  measured  with  the  same  optical  wavefront  instruments  without  distorting  the 
relevant  structures  by  simply  heating  the  wall.  One  consequence  of  this  result  is  the  possibility  of 
using  non-intrusive  optical  diagnostic  tools,  like  a  Malley  probe  or  a  high-speed  wavefront  sensor, 
to  study  boundary  layers  in  low  speed  turbulent  boundary  layer  facilities.  This  result  greatly 
increases  the  number  of  turbulent  boundary  layer  facilities  that  can  be  used  in  aero-optic  research, 
as  well  as  allows  sensitive  optical  sensors  to  be  used  to  study  turbulent  boundary  layers  at  low 
subsonic  speeds;  this  has  already  led  to  a  development  of  new  experimental  approaches  to  study 
fundamental  physics  of  boundary  layers,  described  in  Chapters  4  and  5  of  this  report. 

Experimental  investigations  of  the  effect  of  partial  wall  cooling  have  also  successfully 
demonstrated  notable  aero-optic  mitigation  capabilities,  with  OPD being  reduced  by  as  much 
as  60%  using  this  flow  control  scheme.  This  maximum  reduction  was  achieved  for  cooling 
applied  over  the  full  wall,  and  experiments  showed  that  OPDrms  monotonically  decreases  with 
increasing  RCOoi,  meaning  that  there  is  no  ‘optimal’  partial  wall  cooling  configuration.  Rather,  the 
greatest  reduction  of  aero-optic  aberrations  comes  from  cooling  the  wall  over  the  full  length  of  the 
TBL.  Results  showed,  however,  that  partial  cooling  can  still  provide  significant  reductions  in 
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OPDrms\  cooling  over  64%  of  the  wall  resulted  in  reductions  between  40  -  55%,  while  cooling 
over  just  32%  of  the  wall  reduced  OPDrms  by  about  20  -  40%.  Partial  cooling  results  were  also 
highly  dependent  on  the  location  of  the  applied  cooling,  with  the  largest  reductions  being 
achieved  for  cooling  just  upstream  of  the  optical  aperture  over  longer  lengths  of  the  wall.  This 
result  is  especially  beneficial  for  the  design  of  any  potential  wall  cooling  aero-optic  mitigation 
system  where  cooling  over  the  full  TBL  development  length  is  not  possible  due  to  space  or 
energy  cost  constraints. 

Wall  cooling  at  supersonic  speeds.  Results  presented  in  Section  3.3  indicate  that  the  full 
wall  cooling  model  can  be  easily  extended  for  supersonic  freestream  velocities.  For  partial 
cooling,  Gran,  et  al.  (1974)  demonstrated  that  thermal  sub-layer  effects  are  present  and  similar  for 
partial  cooling  in  supersonic  TBLs.  Therefore,  the  partial  wall  cooling  aero-optic  model  is  also 
expected  to  extend  to  the  supersonic  regime.  One  limitation  of  cooling  for  supersonic  TBL  aero- 
optic  flow  control,  however,  is  found  by  calculating  the  value  of  the  optimal  cooling  temperature. 

As  it  was  shown  in  this  work,  for  full  wall  cooling,  T0ptjmal  =  Taw  - QATJ41 .  Figure  3.11 

presents  estimates  of  T0ptimai  or  several  different  freestream  air  temperatures  that  were  estimated 
from  the  1976  U.  S.  Standard  Atmosphere.  Note  that  in  all  cases  T 'optimal  becomes  equal  to 
absolute  zero  (0  K)  at  about  =  2.13,  which  sets  the  theoretical  maximum  speed  at  which  the 
optimal  wall  cooling  temperature  is  reachable.  Since  cooling  the  wall  to  0  K  is  not  feasible  in 
practical  applications,  however,  it  is  more  practical  to  consider  the  temperatures  of  common 
cryogenic  materials,  such  as  liquid  Oxygen  (LO2),  liquid  Nitrogen  (LN2),  and  liquid  Hydrogen 
(LH2).  Comparing  T0ptimai  to  these  more  realistic  estimates  of  the  minimum  achievable  wall 
cooling  temperature  shows  that  for  LH2  optimal  cooling  temperatures  could  be  achieved  for  < 
2,  and  for  L02  and  LN2  this  range  of  Mach  numbers  shrinks  to  approximately  M, x  <  1.6. 

This  does  not  mean,  however,  that  wall  cooling  would  not  be  an  effective  means  of  aero- 
optic  mitigation  for  supersonic  TBLs  with  M0 0  >  2.  In  fact,  the  aero-optic  model  predicts  that  even 
cooling  the  wall  to  temperatures  that  are  greater  than  Toptimai  will  result  in  considerable  reductions 
of  aero-optic  aberrations.  Estimates  of  the  minimum  value  of  OPDrms  that  can  be  achieved  for 
wall  cooling  temperatures  greater  than  a  set  minimum  value  (i.e.  Tw  >  Toptimai)  are  presented  in 
Figure  3.12  for  LO2,  LN2,  and  LH2.The  results  show  that  even  up  to  Mach  5,  optical  aberrations 
can  be  reduced  by  between  40  -  45%  for  full  wall  cooling  to  temperatures  on  the  order  of  20  - 
100  K.  It  is  also  expected  that  for  partial  wall  cooling,  the  amount  of  aero-optic  be  further 
reduced  as  Rcooi  is  decreased,  as  was  shown  for  subsonic  TBLs.  While  this  is  somewhat  short  of 
the  roughly  60%  reduction  that  was  measured  for  the  subsonic  TBL,  it  is  still  a  significant 
improvement.  Therefore  wall  cooling  should  still  be  considered  as  a  good  candidate  for  use  as  an 
aero-optic  mitigation  technique  in  supersonic  TBLs. 
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Figure  3.11.  Optimal  wall  cooling  temperatures  as  a  function  of  Mach  number  using 
experimentally-obtained  model  constants,  for  different  altitudes  (1976  U.  S.  Standard 
Atmosphere). 


Figure  3.12.  Estimates  of  the  minimum  value  of  OPDrms  that  can  be  achieved  when  Tw  is 
constrained  so  that  Tw  >  Toptimai',  for  full  wall  cooling  with  Tx  =  229  K  (conditions  at 
approximately  30,000  ft.) 
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4.  Aero-Optics  of  Boundary  Layers  at  Low  Reynolds  Numbers 

4.1  Motivation 

As  it  was  discussed  in  the  previous  chapter,  wall  heating  acts  as  a  passive  amplifier  of 
wavefront  distortions.  This  passive  amplification  effect  can  be  used  advantageously  to  study  low 
Reynolds  number  TBLs  that  have  a  weak  aero-optic  signature,  such  as  those  with  low  subsonic 
freestream  velocities  and/or  TBLs  that  are  very  thin.  By  heating  the  wall  it  is  possible  to  a  passive 
‘marker’  to  make  low-Reynolds  number  TBLs  aero-optically  ‘visible.’  Thus  wavefront  sensors 
can  be  used  as  non-intrusive  diagnostic  tools  to  allow  for  investigation  of  the  dynamics  of 
turbulent  boundary  layers  at  relatively  low  Reynolds  numbers.  The  heating  technique  can  also 
provide  valuable  experimental  wavefront  information  that  can  be  used  to  verify/correct  various 
scaling  models  at  low  Reynolds  numbers  and  to  provide  a  direct  comparison  with  computational 
simulations,  as  currently  there  is  a  large  Reynolds  number  ‘gap’  between  experimental 
measurements  with  typical  Ree>  20,000  (Gordeyev  et  al,  2014)  and  numerical  simulations  (Wang 
&  Wang,  2012;  White  &  Visbal,  2012)  with  much  lower  Re#~  1,000  to  3,000. 

4.2  Results 

One  set  of  measurements  were  conducted  in  the  Hessert  Transonic  Wind  Tunnel  (TWT)  at 
the  University  of  Notre  Dame.  The  total  length  to  the  beginning  of  the  optical  section  was 
shortened  to  approximately  100  cm  to  reduce  the  boundary  layer  thickness,  and  therefore  the 
Reynolds  number.  Non-adiabatic  wall  temperatures  were  introduced  by  replacing  the  Plexiglas 
wall  with  an  aluminum  wall  with  an  8  mm  thick  Aluminum  plate  for  the  first  100  cm  of  the  test 
section,  and  heating  it  from  the  outside  surface  using  strips  of  flexible  electric  resistive  coil 
heaters.  Just  downstream  of  the  boundary  layer  development/wall  heating  section,  portions  of  the 
Plexiglas  on  the  upper  and  lower  walls  of  the  wind  tunnel  were  replaced  with  optical  quality  glass 
plates  to  allow  accurate  optical  characterization  of  the  boundary  layer  using  aero-optic  wavefront 
sensors.  A  Pitot  tube  mounted  downstream  of  the  optical  measurement  location  was  used  to 
measure  the  ffee-stream  velocity  throughout  the  experiment. 

Additional  experiments  were  performed  at  California  Institute  of  Technology  in  the 
Merrill  Wind  tunnel  (MWT),  a  recirculating  tunnel  with  a  0.6x0. 6x2.4  m  constant  area  test 
section.  The  wall  of  a  test  plate  mounted  in  the  center  of  the  test  section  was  heated  at  two 
streamwise  locations  using  the  same  flexible  heating  panels  as  in  the  Notre  Dame  experiments, 
this  time  sandwiched  between  the  aluminum  test  surface  and  a  heavily  insulated  lower  surface. 


Table  4.1  Boundary  Layer  Wall  Heating  Conditions 


Facility 

Foo  [m/s]  M 

S  [cm] 

Re# 

\T  [K] 

./samp  [kHz] 

Caltech  MWT 

9.4 

0.03 

2.7 

1,700 

21 

30 

ND  TWT 

64.8 

0.18 

1.2 

4,200 

15-28 

200 

11 

98.8 

0.28 

1.2 

5,700 

9-24 

200 

11 

118 

0.35 

1.2 

7,900 

7-28 

200 

11 

140 

0.41 

1.2 

9,000 

7-24 

200 

11 

140 

0.41 

2.4 

20,000 

0 

200 

Wavefront  measurements  were  acquired  downstream  of  the  heated  wall  boundary  layer 
development  section  in  both  facilities  using  the  Malley  Probe  1-D  wavefront  sensor.  The  Mach 
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number  M,  boundary  layer  thickness  d,  Reynolds  number  Re#,  wall  temperature  A  T,  and  sampling 
rate  /sam P,  corresponding  to  each  of  these  cases  are  given  in  Table  4.1.  From  the  deflection  angle 
time  series  measured  with  the  Malley  probe,  deflection  angle  amplitude  spectra  and  OPDrms  were 
computed. 

Figure  4.1  presents  normalized  deflection  angle  spectra  obtained  using  wall  heating  as  a 
passive  amplifier  for  the  different  Reynolds  numbers  spanning  an  order  of  magnitude,  from  Reg  = 
1,700  to  20,000.  In  the  vicinity  of  the  spectra  peak,  and  for  the  low  end  of  the  spectral  peak  ( Sts< 
1),  all  cases  show  good  collapse.  The  peaks  for  all  measured  spectra  are  at  approximately  St\j  = 
1.0,  independent  of  Reynolds  number.  In  Gordeyev,  et  al.  (2014)  it  was  shown  that  using 
Kolmogorov-like  arguments  for  the  inertial  range,  the  deflection  angle  amplitude  spectrum  at 
large  frequencies  should  behave  as  ~/(_2/3) ;  this  curve  is  plotted  in  Figure  4.1  as  a  dashed  black 
line.  The  spectrum  at  the  lowest  Reynolds  number  quickly  falls  off  of  this  theoretically  predicted 
behavior,  indicating  a  fairly  small  inertial  range,  while  the  fall-off  is  less  drastic  for  larger 
Reynolds  numbers,  implying  a  larger  inertial  range. 


Figure  4.1.  Normalized  deflection  angle  spectra  for  different  Reynolds  numbers. 


To  model  the  observed  changes  in  the  (high  frequency  slope)  size  of  the  inertial  range  of 
deflection  angle  spectra  for  different  Reynolds  numbers,  the  spectral  model  from  presented  in 
Gordeyev,  et  al.  (2014)  was  modified  in  order  to  account  for  changes  in  the  spectra  roll-off  as  a 
function  of  Re,?.  The  fall-off  of  the  spectra  was  assumed  to  take  a  form  inspired  by  Tatarski’s 
modification  of  Kolmogorov’s  atmospheric  wavefront  spectrum  to  account  for  the  presence  of 
inner  scale  dissipative  structures  (Tatarski,  1971).  Seeking  a  function  in  the  form  of  the  original 
model  times  the  exponential  term,  the  modified  spectral  model  was  found  to  be 

r  o.  \2' 


St  c 


peak 


l  +  (SG/0.83)5/3 


exp 


St, 


/(Re„). 


(4.1) 


The / (Re,?)  term  was  empirically  determined  to  be / (Re,?)  ~  1.6Re^(0  22)  by  fitting  equation  (4.1)  to 
the  experimental  data  presented  in  Figure  4.1. 
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Ree 

Figure  4.2.  Empirically  determined  values  of /  (Re#)  from  equation  (4.1)  and  the  corresponding 
power  law  fit. 


Figure  4.3  shows  comparisons  between  experimental  data  and  the  modified  model, 
equation  (4.1),  for  four  Reynolds  numbers  ranging  from  1,700  to  20,000.  For  Stg<  0.5,  the 
experimentally  measured  spectra  and  model  spectra  deviate  from  one  another,  with  the 
experimentally  measured  spectra  containing  more  energy.  This  behavior  is  consistent  with  results 
from  previous  studies,  where  it  was  shown  that  for  Malley  probe  measurements,  this  increase  in 
energy  is  a  result  of  contamination  from  mechanical  vibrations  in  Malley  probe  data  that  cannot 
be  properly  removed  (Gordeyev  et  al,  2014).  For  Stg>  0.5,  Figure  4.3  shows  that  the  modified 
spectral  model,  equation  (4.1),  does  a  good  job  of  describing  the  high-frequency  roll-off  behavior 
for  a  large  range  of  Reynolds  numbers. 

Heated  wall  TBL  wavefront  statistics  can  also  be  used  to  evaluate  the  validity  of  scaling 
laws  for  OPDrms  at  low  Reynolds  numbers.  One  method  of  doing  this  is  to  rearrange  the  terms  of 
the  TBL  aero-optic  model  presented  in  equation  (3.18)  to  solve  for  the  model  constant,  B: 


B  =  OPDrj£T  =  Q) 
KGDpJM2^Cf 


(4.2) 


Here  the  only  term  that  explicitly  depends  on  Reynolds  number  is  C/,  and  there  are  a  number  of 
relations  for  computing  C/  as  a  function  of  Re#  that  compare  with  experimental  measurements 
(Nagib,  et  al.,  2007).  Therefore  if  the  assumptions  used  to  derive  equation  (3.18)  are  valid  at  low 
Reynolds  numbers,  the  model  constant  should  be  consistent  with  what  was  found  at  high  Re#, 
where  B  =  0.19  ±  0.01  (Gordeyev  et  al.,  2014). 

Equation  (4.2)  was  used  to  estimate  B  from  experimental  data  obtained  at  low  Reynolds 
number,  using  values  of  OPDrms  (AT  =  0)  measured  directly  for  adiabatic  wall  conditions  and 
values  of  OPDrms  (AT  =  0)  estimated  from  a  linear  fit  of  heated  wall  data.  The  results  of  these 
estimates  are  shown  in  Figure  4.4,  compared  with  the  previously  obtained  experimental  value  of 
B  =  0.19±0.01  from  Gordeyev,  et  al.  (2012).  Good  agreement  was  found  between  values  of  B 
estimated  from  equation  (4.2)  for  both  adiabatic  wall  wavefront  measurements  and  heated  wall 
wavefront  measurements  at  low  Reynolds  numbers,  especially  for  Re  >  5,700.At  Re  =  4,200, 
there  is  a  slight  increase  in  the  value  of  Bestimated  to  about  0.22  -  0.23,  however  the  error  bars  for 
these  estimates  are  large. 
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St<5  Stj 

Figure  4.3.  Comparison  of  experimental  deflection  angle  spectra  (blue  line)  and  the  modified 
empirical  model  from  equation  (4.1)  (red  line)  for  Reynolds  number,  Re<?,  a)  1,700  (Caltech),  b) 
4,200  (Notre  Dame),  c)  5,700  (Notre  Dame),  and  d)  20,000  (Notre  Dame). 

4.3  Conclusions  and  Discussion 

This  chapter  presented  the  results  of  aero-optic  wavefront  measurements  of  heated-wall 
boundary  layers  for  a  range  of  subsonic  Mach  numbers  from  0.03  to  0.4,  and  Reynolds  numbers 
Re#  =  1,700  -  20,000.  The  effect  of  increased  wall  temperatures  on  deflection  angle  spectra  was 
analyzed,  and  it  was  found  that  heating  indeed  acted  as  a  passive  amplifier  on  TBL  turbulence, 
confirming  that  the  wall  heating  scaling  laws  developed  by  Cress  (2010)  for  higher  Reynolds 
number  TBLs  were  valid  for  Reynolds  numbers  at  least  as  low  as  4,200.  Using  the  scaling  laws,  it 
was  shown  that  equation  (3.18)  for  heated  walls  could  be  used  to  re-scale  heated- wall  wavefront 
data  to  recover  wavefront  statistics  for  adiabatic  wall  temperatures  at  low  Reynolds  numbers. 
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Figure  4.4.  Model  constant  B  from  equation  (3.18)  estimated  from  low-Reynolds  number 
wavefront  measurements  from  OPDrms. 


Wavefront  spectra  from  both  heated  and  un-heated  walls  were  extracted  and  compared  to 
investigate  the  effect  of  Reynolds  number  on  the  shape  of  boundary  layer  wavefront  spectra.  It 
was  found  that  the  inertial  sub-range  of  aero-optic  wavefront  spectra  is  reduced  with  decreasing 
Reynolds  number,  which  is  consistent  with  previous  observations  of  other  turbulence  spectra 
(Tatarski,  1971).  The  empirical  spectral  model  was  modified  in  order  to  model  this  Reynolds 
number  dependence  by  the  addition  of  an  exponential  decay  term  (which  is  a  function  of 
Reynolds  number)  that  approximates  the  high  frequency  roll-off  observed  in  the  experimentally 
measured  wavefront  spectra.  The  modified  model  was  shown  to  agree  well  with  spectra  obtained 
from  experiments  spanning  an  order  of  magnitude  in  Reynolds  number,  from  Re#  =  1,700  - 
20,000. 

The  presented  method  for  extracting  wavefront  statistics  TBLs  from  heated  wall 
wavefront  data  provides  a  novel  and  promising  method  for  studying  the  aero-optic  characteristics 
of  TBLs  that  would  otherwise  be  difficult  to  measure  due  to  the  problem  of  aero-optic 
‘invisibility’  that  occurs  when  Mach  and/or  Reynolds  numbers  in  the  TBL  are  low.  This  approach 
for  obtaining  low  Reynolds/Mach  number  TBL  wavefront  data  is  especially  useful  for  studying 
low-Reynolds  number  TBLs  using  non-intrusive  aero-optic  wavefront  sensors.  This  technique, 
along  with  the  modification  proposed  to  the  spectral  model  for  Reynolds  number  dependence,  are 
important  developments  in  obtaining  experimentally  measured  low  Reynolds  number  TBL 
wavefront  statistics  for  comparison  to  computational  results,  which  are  limited  to  low  Reynolds 
numbers  by  limitations  of  practical  parameters  for  simulations. 
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5.  Simultaneous  Optical/Velocity  Measurements 


5.1  Motivation 

Large-scale  structures,  located  in  the  outer  part  of  the  subsonic  boundary  layer  carry  an 
important  role  in  boundary  layer  dynamics,  responsible  for  entrainment  process  and,  via  a  link 
with  small-scale  structures  near  the  wall,  for  instantaneous  drag  near  the  wall,  to  mention  a  few 
(Marusic  et  al.,  2010).  While  there  exists  a  large  body  of  experimental  research  about  the  large- 
scale  structures  (Smits  et  al,  2011),  there  are  open  questions  about  their  topology,  dynamics  and 
the  interaction  with  the  near-wall  structures. 

Traditionally,  large-scale  structures  are  characterized  by  the  velocity  field,  which  is 
measured  using  hot-wires  or  a  PIV  technique.  These  techniques  give  either  detailed  temporal 
information  in  a  few  spatial  points  (hot-wires)  or  potentially  time-resolved  spatial  velocity  field 
information  in  a  plane  (particle  image  velocimetry,  PIV).  In  order  to  get  instantaneous  three- 
dimensional  information  about  the  structure,  one  has  to  resort  to  rather  complex  techniques,  like 
tomographic-PIV  (Elsinga  et  al.,  2006)  or  plenoptic  (Farbinger  and  Thurow,  2012)  techniques,  for 
example.  An  alternative  way  to  non-intrusively  study  the  characteristics  of  large-scale  structures 
in  boundary  layers  is  to  measure  related  density  distortions  uisng  various  wavefront  sensors. 

5.2  Approach 

In  this  work,  we  leverage  both  velocimetry  and  optical  techniques  to  characterize  large 
scale  structure  in  incompressible  and  compressible  turbulent  boundary  layers. . 

As  stated  before,  in  incompressible,  wall-heated  boundary  layers  density  fluctuations 
occur  due  to  total  temperature  variations.  Above  M  =  0.3,  compressibility  effects  can  also  change 
the  density.  If  the  total  temperature  is  assumed  to  be  a  constant,  which  is  equivalent  to  zero 
pressure  fluctuations,  the  density  fluctuations,  p\  in  boundary  layers  are  due  to  the  static 
temperature  fluctuations,  T,  (adiabatic  cooling/heating),  which  is  the  base  for  the  Strong 
Reynolds  Analogy,  discussed  earlier  in  Section  3.2. 

The  time-averaged  version  of  SRA  was  shown  to  correctly  estimate  the  time-averaged 
levels  of  BL  aero-optical  distortions  (Gordeyev  et  al.  2014).  If  instantaneous  wavefronts  and 
velocity  fields  are  known,  it  is  possible  to  estimate  the  validity  of  the  instantaneous  version  of 
SRA, 

CpT'  (x,  t)  =  —U u'  (x,  t)  (5.1) 

by  comparing  measured  wavefronts  and  estimated  wavefronts  from  instantaneous  velocity  field 
using  equations  (1.1)  and  (5.1). 

In  Chapter  3  it  was  shown  that  if  the  boundary  layer  wall  is  moderately-heated  relative  to 
the  freestream  density,  it  will  not  modify  the  underlying  velocity  structure.  Instead,  moderate 
heating  will  simply  introduce  passive  temperature  markers  in  the  boundary  layer  and  effectively 
amplify  aero-optical  levels.  In  other  words,  by  introducing  moderate-temperature  mismatch 
between  the  wall  and  the  boundary  layer,  it  is  possible  to  thermally  tag  and  measure  three- 
dimensional  large-scale  structures  in  incompressible  boundary  layers. 

The  similarities  between  compressible  and  incompressible  aero-optical  distortion  statistics 
have  two  implications.  First,  they  suggest  that  aero-optical  distortion  studies  can  be  carried  out  in 
incompressible  tunnels,  opening  up  more  avenues  for  measurement.  Secondly,  and  arguably  more 
importantly,  they  suggest  that  inherent  features  of  turbulence,  rather  than  specific  compressibility 
effects,  underlie  the  patterns  of  aero-optical  distortion  seen  in  compressible  and  heated, 


60 


incompressible  flow.  With  these  ideas  in  mind,  this  study  examines  aero-optical  distortion  in  a 
heated,  incompressible  flow  while  simultaneously  measuring  flow  features  to  examine  how 
turbulent  structures  affect  laser  distortion. 

Two  simultaneous  wavefront-velocity  experiments  were  conducted  at  a  low  subsonic 
speed  at  California  Institute  of  Technology  (Caltech),  using  the  moderately-heated  approach  and 
at  a  compressible  subsonic  speed  at  University  of  Notre  Dame.  A  Malley  probe  was  used  in  the 
Caltech  experiments  to  characterize  optical  distortions  and  obtain  an  appropriate  mean  convection 
velocity.  A  Shack-Hartmann  sensor  was  used  to  obtain  more  detailed,  two-dimensional  aero- 
optical  information  in  the  Notre  Dame  experiments,  specifically,  the  two-dimensional,  time- 
resolved  OPD. 

5.3  Experimental  Setup:  Caltech 

Experiments  are  conducted  in  the  Merrill  wind  tunnel  at  Caltech,  a  recirculating  facility 
with  a  2ft  x  2ft  test  section.  A  turbulent  boundary  layer  is  developed  over  a  partially  heated  flat 
plate  at  a  Reynolds  number  of  3300.  The  plate  is  heated  to  between  15  and  20  C  above  the  free 
stream  temperature.  A  Malley  probe  is  used  to  measure  aero-optical  distortion  and  Particle  image 
velocimetry  (PIV)  is  used  to  measure  the  velocity  field  in  a  wall-normal,  streamwise  plane.  The 
PIV  laser  sheet  is  aligned  to  be  coincident  with  the  Malley  probe  beams  in  the  spanwise  and 
streamwise  directions.  PIV  images  are  taken  at  1500  Hz,  double  pulsed  with  a  35  ps  separation 
between  pulses.  The  field  of  view  is  60  mm  x  60  mm  (1.7c>  x  \.15).  A  schematic  of  the 
experimental  setup  showing  the  thermal  boundary  layer  the  passage  of  the  Malley  probe  beams 
through  the  flow,  and  the  location  of  the  PIV  light  sheet  in  the  flow  is  presented  in  Figure  5.1. 


Figure  5.1.  Schematic  of  the  experimental  setup.  A  turbulent  boundary  layer  is  shown  developing 
over  a  partially  heated  flat  plate.  A  thermal  boundary  layer  develops  with  an  inner  and  outer  edge. 
The  Malley  probe  beams  are  passed  up  and  back  through  the  flow.  The  red  line  indicates  the 
upstream  beam  while  the  blue  indicates  the  downstream  beam.  The  green  triangle  indicates  the 
location  of  the  PIV  light  sheet  in  the  flow. 

The  heated  plate  that  seeds  density  variation  into  the  flow  ends  upstream  of  the 
measurement  location  creating  an  inner  thermal  boundary  layer.  The  temperature  was  measured 
with  a  cold  wire  as  a  function  of  height  over  a  fully  heated  wall  and  over  the  measurement 
location.  Figure  5.2  shows  these  temperature  profiles  and  compares  them  to  data  from  Antonia& 
Kawamura,  1977  of  a  heated  turbulent  boundary  layer.  The  fully  heated  and  partially  heated 
profiles  are  seen  to  diverge  at  a  y+  =  100  or  y 18=  0.13.  The  partial  heating  allows  for  a  focus  in 
the  outer  boundary  layer,  which  highlights  the  regions  of  interest  in  the  flow  for  this  study.  It  also 
does  not  seem  to  impact  the  comparison  to  the  compressible  case,  as  the  average  convection 
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velocities  and  spectra  look  the  same  between  the  partially  heated,  incompressible  case  and  the 
compressible  case. 


Figure  5.2.  Normalized  averaged  temperature  relative  to  the  wall  temperature  as  a  function  of 
distance  from  the  wall.  The  temperature  profiles  that  plateau  are  measured  at  the  measurement 
location  (A)  while  the  temperature  profiles  that  decay  to  zero  are  measured  directly  over  the 
heated  wall  (B).  The  average  temperature  profile  in  the  partially  heated  case  diverges  from  the 
fully  heated  case  at  roughly  y+  =  100  or  y/S=  0.13  (Rought  2013;  Antonia  et  al  1977). 


Sample  data  from  the  Malley  probe  and  PIV  measurements  are  shown  in  Figure  5.3.  The 
streamwise  deflection  angles  of  the  two  Malley  probe  beams  are  shown  as  a  function  of  time  in 
Figure  5.3(a).  The  vertical  black  line  indicates  the  time  at  which  the  PIV  data  shown  in  Figure 
5.3(b)  was  taken.  Figure  5.3(b)  shows  a  vector  plot  of  the  velocity  field  with  a  convection  velocity 
of  0.83Uoo  subtracted  off  to  show  the  structure  in  the  flow.  The  locations  of  the  two  Malley  probe 
beams  are  indicated  by  the  red  (upstream  beam)  and  blue  (downstream  beam)  lines.  The  figure 
shows  a  smaller  field  of  view  than  is  available  in  order  to  more  easily  see  the  structure. 

5.4  Results:  Caltech 

5.4.1  Statistics 

The  mean  and  RMS  velocity  statistics  are  found  from  the  PIV  data  and  compared  to  data 
from  DeGraaff  &  Eaton,  2000  at  similar  Reynolds  numbers.  The  data  is  found  to  agree  well  above 
a  y+  of  about  40  or  equivalently  y 18=  0.05.  Figures  5.4(a)  and  5.4(b)  show  the  mean  and  RMS 
compared  to  the  DeGraaff  &  Eaton,  2000  data. 
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Figure  5.3.  Sample  Malley  probe  (a)  and  PIV  (b)  data.  Red  line  indicates  upstream  Malley  probe 
beam  while  blue  line  indicates  downstream  Malley  probe  beam,  (a)  shows  the  streamwise  angle 
as  a  function  of  time  for  each  beam.  The  black  vertical  line  in  (a)  indicates  the  time  at  which  the 
PIV  data  shown  in  (b)  was  taken,  (b)  shows  the  instantaneous  velocity  field  with  0.83  Uoo 
subtracted  off  to  be  able  to  see  structure.  The  field  of  view  shown  is  smaller  than  the  available 
field  of  view  to  more  easily  see  structure  in  the  flow.  Vertical  red  and  blue  lines  are  drawn  in  to 
indicate  the  locations  of  the  upstream  and  downstream  Malley  probe  beam  paths. 
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Figure  5.4.  Mean  profile  (a)  and  RMS  (b)  as  a  function  of  distance  from  the  wall.  Blue  points 
indicate  data  from  PIV  measurements  in  this  study.  Black  and  red  points  indicate  data  from 
DeGraaff  &  Eaton,  2000  at  similar  Reynolds  numbers  for  comparison.  Data  is  seen  to  agree  above 
a  y+  of  about  40  and  agree  very  well  above  a  y+  of  about  100  (DeGraaff  &  Eaton,  2000). 

5.4.2  Convection  Velocities 

Convection  velocity  information  is  obtained  from  the  aero-optical  distortion  measurement 
using  a  Pearson  correlation,  defined  as: 
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where  0i  is  the  upstream  beam  and  02  is  the  downstream  beam.  The  value  of  x  that  gives  the 
largest  value  of  R  is  taken  as  the  convection  time  for  a  disturbance  to  travel  from  the  upstream 
Malley  probe  beam  to  the  downstream  Malley  probe  beam. 

The  correlation  can  be  done  over  small  periods  of  time  (in  this  study,  local  correlation  is 
performed  over  1.4  ms  or  0.5  convective  time  units)  to  find  a  local  convection  velocity  that 
ideally  corresponds  to  a  single  structure,  or  can  be  done  over  the  full  length  of  the  signal  to  find  a 
global  convection  velocity  of  the  structures.  The  average  of  all  local  convection  velocities  is 
found  to  be  0.85Uoo  with  a  standard  deviation  of  O.lUoo.  The  global  convection  velocity  is  found  to 
be  0.83Uoo.  These  convection  velocities  of  about  0.8  U®  suggest  that  something  in  the  outer 
boundary  layer  is  most  responsible  for  the  deviations  in  the  laser  path.  By  comparing  the 
convection  velocity  to  the  mean  velocity  profile,  an  average  height  of  the  relevant  structures  can 
be  computed.  Figure  5.5  shows  the  heights  corresponding  to  0.83Uoo  and  0.85Uoo  in  an 
instantaneous  velocity  field  (a)  and  the  mean  profile  (b).  An  average  height  of  OAU®  is  observed. 
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Figure  5.5.  The  height  corresponding  to  the  average  convection  velocity  of  angle  deviations  is 
indicated  relative  to  an  instantaneous  vector  plot  (a)  and  the  mean  velocity  profile  (b).  The  solid 
green  line  indicates  the  height  at  which  the  mean  velocity  is  0.85  Uoo,  the  average  value  of  the 
local  convection  velocities.  The  dashed  green  line  indicates  the  height  of  at  which  the  mean 
velocity  is  0.83  Uoo,  the  global  convection  velocity  found  by  correlating  the  full  signal. 


5.4.3  Conditional  Averaging 

Conditional  averaging  is  next  used  to  find  patterns  in  the  velocity  field  that  consistently 
leads  to  large  deflections  in  the  laser  path.  The  velocity  field  surrounding  a  Malley  probe  laser 
path  is  conditionally  averaged  on  large  positive  (defined  as  downstream)  deflections  and  large 
negative  (upstream)  deflections  separately.  Figure  5.6(a)  shows  the  average  velocity  field  given  a 
condition  of  a  positive  deviation  in  the  laser  angle  larger  than  half  of  the  maximum  deviation, 
while  Figure  5.6(b)  shows  the  average  velocity  field  for  a  negative  deviation  in  the  laser  angle 
less  than  a  third  of  the  minimum  deviation.  The  condition  on  positive  deviations  is  stronger  than 
that  on  negative  in  order  to  have  roughly  50  frames  averaged  per  1024  frame  run  for  both  positive 
and  negative  averages.  The  difference  in  the  conditions  is  caused  by  a  non-uniform  distribution  in 
laser  deflections,  reflected  in  the  negative  skewness  of  the  deflection  angle  signal. 

Each  figure  has  a  particular  convection  velocity  subtracted  from  the  full  flow,  which 
highlights  a  particular  region  of  the  flow.  The  average  of  the  local  convection  velocities  found  for 
the  large  positive  deflections  is  subtracted  from  the  full  velocity  field  in  Figure  5.7.  The 
equivalent  procedure  is  done  for  the  large  negative  deflection  convection  velocities  for  Figure  5.8. 
For  large  positive  excursions  in  the  laser  deflection  angle,  a  vortex  pattern  is  observed  in  the 
conditionally  averaged  flow,  while  for  large  negative  excursions,  a  stagnation  point  is  observed. 
Both  are  observed  to  travel  with  the  same  approximate  convection  velocity  of  0.8  5  Uoo  and  at  a 
height  0.55. 

However,  if  the  subtracted  convection  velocity  is  changed,  the  same  feature  is  visible  at 
different  heights  in  the  flow.  Figures  5.7  and  5.8  show  the  change  in  position  of  the  vortex  and 
stagnation  point  respectively  as  the  subtracted  convection  velocity  is  changed.  The  conditionally 
averaged  flow  highlights  a  feature  that  occurs  at  many  heights  and  convection  velocities  in  the 
flow.  This  suggests  that  the  laser  deflection  angle  can  be  affected  by  structures  in  a  relatively 
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large  region  of  the  flow,  but  the  relationship  between  the  laser  deflection  and  shape  of  the 
structure  is  consistent.  Vorticity  along  the  beam  path  generates  downstream  deflection,  while 
stagnation  and  shear  along  the  beam  path  generates  upstream  deflection. 


(a)  (b) 

Figure  5.6.  Average  of  the  velocity  field  around  the  laser  beam  path  during  large  excursions  in  the 
laser  deflection  angle,  0.  Large  positive  excursions,  defined  as  0  >  0.5  max  0,  produce  a  vortex 
structure  in  the  velocity  field  (a).  Large  negative  excursions,  defined  as  0  <  l/3*min  0,  produce  a 
stagnation  point  in  the  velocity  field  (b).  Positive  0  is  defined  as  a  downstream  bend  in  the  laser, 
while  negative  0  is  defined  as  an  upstream  bend.  The  average  convection  velocity  of  the 
deflection  events  is  subtracted  from  the  velocity  field,  in  this  case  0.85  Uoo,  allowing  the  structures 
to  be  visible.  The  red  line  indicates  the  position  of  the  Malley  probe  beam  in  the  velocity  field. 
Both  structures  are  observed  to  sit  at  a  height  of  0.55. 
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Figure  5.7.  Vector  plots  show  a  close-in  perspective  on  the  conditionally  averaged  velocity  data 
shown  in  Figure  5.6a  with  varying  convection  velocities  subtracted  off.  Each  velocity  allows  a 
vortex  to  be  visible  at  a  different  height  in  the  flow.  The  presence  of  vorticity  along  the  beam  path 
in  a  large  region  of  the  flow  is  consistent  with  strong  downstream  deflections  in  the  beam  path. 
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Figure  5.8.  Vector  plots  show  a  close-in  perspective  on  the  conditionally  averaged  velocity  data 
shown  in  Figure  5.6b  with  varying  convection  velocities  subtracted  off.  Each  velocity  allows  a 
stagnation  point  to  be  visible  at  a  different  height  in  the  flow.  The  presence  of  a  stagnation  point 
along  the  beam  path  in  a  large  region  of  the  flow  is  consistent  with  strong  upstream  deflections  in 
the  beam  path. 

5.4.4  Correlation 

To  continue  probing  the  relationship  between  the  velocity  field  and  the  laser  distortion, 
correlations  are  carried  out  between  the  laser  probe  deflection  angle  and  the  velocity  at  a  specific 
point  in  space  for  each  PIV  frame.  The  vertical  velocity  is  extracted  from  velocity  field  at  the  x- 
location  of  the  Malley  probe  beam  and  a  y-location  yc.  yc  is  chosen  using  convection  velocity 
arguments  demonstrated  in  Figure  5.5.  The  local  convection  velocity  of  the  Malley  probe  signal 
is  compared  to  the  mean  velocity  profile.  The  height  at  which  the  mean  profile  matches  the  local 
convection  velocity  is  taken  as  yc. 

The  vertical  velocity,  v(xMp,yc),  is  correlated  to  the  laser  deflection  angle  x  at  each  PIV 
time  step.  A  Pearson  correlation  coefficient,  defined  as: 
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of  R  =  -0.28  is  found  with  a  time  delay  normalized  by  outer  units  of  xLE/S  =  0.12.  Figure  5.9 
shows  the  correlation  coefficient  as  a  function  of  the  time  delay  x. 
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Figure  5.9.  Correlation  of  the  vertical  velocity  v  with  0  as  a  function  of  time  delay.  The  maximum 
absolute  value  of  the  correlation  coefficient  is  found  at  a  normalized  time  delay  of  0.13.  The 
negative  sign  of  the  correlation  indicates  that  opposite-signed  deviations  in  v  and  6  are  related. 
The  time  delay  suggests  that  the  part  of  the  feature  that  has  a  strong  velocity  disturbance  and  the 
part  that  strongly  affects  the  laser  propagation  angle  are  separated  with  some  specific  spacing  in 
the  flow. 


(a)  (b) 

Figure  5.10.  A  comparison  of  conditional  averaging  and  correlation  results.  The  vector  plots 
shown  are  the  conditional  averages  on  large  positive  (a)  and  negative  (b)  deviations  in  0.  The  red 
line  indicates  the  location  of  the  0  event.  The  black  dashed  line  shows  the  x-location  where  there 
should  be  a  correlated  negative  (a)  or  positive  (b)  v  event,  given  the  time  delay  found  in  the 
correlation.  At  a  height  of  0.55  the  velocity  fields  show  a  strong  -v  event  (a)  and  strong  +v  event 
(b)  at  the  x-location  where  the  correlation  predicted  these  events,  indicating  that  the  conditional 
averaging  and  correlation  results  agree. 

The  negative  correlation  indicates  that  opposite-signed  deviations  in  v  and  6  are  related. 
The  time  delay  suggests  that  the  part  of  the  flow  that  has  a  strong  velocity  disturbance  occurs 
upstream  of  the  part  of  the  flow  that  causes  a  strong  deflection  in  the  laser  path.  In  order  to 
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investigate  this  aspect  of  the  correlation,  the  time  delay  is  translated  to  a  spatial  displacement 
using  the  local  convection  velocity.  That  spatial  displacement  is  then  compared  to  the  results  of 
the  conditional  averaging.  The  location  of  correlated  events  is  overlaid  with  the  conditional 
averaging  results  in  Figures  5.10(a)  and  5.10(b).  The  vector  field  in  each  plot  is  the  conditional 
average  shown  previously.  The  red  line  indicates  the  location  of  a  large  deviation  in  laser  angle, 
while  the  black  line  indicates  the  location  of  a  correlated  velocity  event  of  the  opposite  sign.  Good 
agreement  is  seen  between  the  location  of  the  black  line  and  strong  velocity  signals  of  the 
appropriate  sign  in  the  conditionally  averaged  velocity  fields. 

5.5  Discussion 

The  presence  of  structure  in  the  conditionally  averaged  velocity  field  and  the  correlation 
between  the  vertical  velocity  and  the  laser  deflection  both  support  the  idea  that  structures  in 
turbulence  shape  the  laser  deflection  angle.  However,  the  mechanism  of  how  these  structures 
create  these  laser  deflections  is  perhaps  not  obvious  from  the  data.  A  hypothesis  is  proposed  that 
fits  the  patterns  of  the  data,  and  opens  room  for  more  questions. 

The  vortex  and  stagnation  point  seen  in  the  conditionally  averaged  velocity  fields  are 
hypothesized  to  be  associated  with  hairpin  vortices.  Hairpins  have  been  studied  at  great  length  by 
the  turbulence  community  and  have  the  features  that  are  seen  in  the  conditionally  averaged  flow 
field.  They  sit  in  the  appropriate  region  of  the  flow  to  have  an  average  convection  velocity  of 
O.8U00  and  the  hairpin  packets  are  of  the  right  size  to  give  a  peak  in  the  deflection  spectrum  at  St§ 
of  order  one  (Adrian  et  al  2000). 

Hairpins  have  previously  been  argued  to  lie  along  an  interface  between  fast  and  slow 
moving  'momentum  zones'  and  form  hairpin  packets  (Adrian  et  al  2000).  The  fast  and  slow 
moving  fluid  may  be  associated  with  cold  (free  stream)  and  warm  (near-wall)  fluid  respectively 
using  passive  scalar  transport  arguments.  If  hairpins  lie  along  a  boundary  between  warm  and  cold 
regions  of  air,  they  might  shape  that  boundary  in  ways  that  consistently  deflect  the  laser,  giving 
the  results  shown.  Figure  5.11  shows  a  model  of  the  separation  between  warm  and  cold  regions  of 
flow.  The  specific  gradient  at  the  vortex  head  that  causes  a  upstream  leaning  deviation  requires 
further  study. 


Figure  5.11.  Model  of  warm  (red)  and  cold  (blue)  regions  of  the  flow.  Grey  circles  represent 
hairpin  heads,  with  the  shear  layers  trailing  down  the  backside  of  the  temperature  interface.  Laser 
paths  are  drawn  in  to  represent  the  types  of  deflections  seen  in  the  conditional  averaging.  The 
slight  gradient  in  color  at  the  bottom  of  the  warm  region  indicates  the  height  cold  inner 
temperature  region  present  due  to  the  wall  heating  element  ending  before  the  measurement 
location. 
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The  observation  that  structures  in  turbulence  are  tied  directly  to  laser  deflections  opens 
many  new  questions.  Future  work  will  focus  on  how  the  structures  shape  the  density  field  and 
therefore  affect  the  laser  deflections,  as  well  as  how  the  structures  may  be  tied  to  the  pickup  of 
heat  at  the  wall.  Additionally,  future  work  will  explore  the  range  of  sizes  and  convection 
velocities  associated  with  optically  aberrating  structures  to  better  quantify  the  scales  that  are 
optically  important  in  turbulent  boundary  layers. 

5.6  Conclusions:  Caltech 

This  study  aims  to  understand  the  connection  between  aero-optical  distortion  in  turbulent 
flows  and  turbulent  velocity  structure.  Aero-optical  distortion  and  the  fluid  velocity  field  are 
simultaneously  measured  in  a  heated  turbulent  boundary  layer.  Conditional  averages  are  used  to 
determine  the  velocity  structure  that  yields  large  deflections  in  the  laser  propagation  angle.  Large 
downstream  deflections  in  the  laser  are  found  to  be  consistent  with  the  presence  of  a  vortex 
passing  through  the  laser  path.  Large  upstream  deflections  in  the  laser  are  found  to  be  consistent 
with  the  presence  of  a  stagnation  point  or  shear  layer  passing  through  the  beam  path.  A  Pearson 
correlation  is  performed  between  the  vertical  velocity,  v,  and  the  laser  deflection  angle,  yielding  a 
correlation  coefficient  of  R  =  -0.28.  A  specific  separation  in  space  between  the  velocity  and 
deflection  angle  events  is  found  to  maximize  their  correlation,  suggesting  a  specific  size  of  the 
turbulent  structure.  A  hypothesis  is  put  forward  that  the  beam  is  distorted  by  the  passage  of 
hairpin  vortices,  which  are  characterized  by  both  vortex  heads  and  trailing  shear  layers.  The 
hairpins  are  proposed  to  lie  along  the  interface  between  warm  and  cool  regions  in  the  flow,  much 
as  they  have  been  shown  to  lie  along  the  interface  between  regions  of  low  and  high  streamwise 
momentum.  Hairpin  vortices  are  thought  to  affect  the  beam  deflection  angle  by  shaping  a  sharp 
gradient  in  density  in  the  flow. 

5.7  Experiments:  Notre  Dame 

Measurements  were  obtained  in  the  Hessert  Transonic  Wind  Tunnel  (TWT)  at  the 
University  of  Notre  Dame.  Freestream  velocity  was  held  constant  at  M  =  0.4  for  the  duration  of 
the  measurements. 


Figure  5.12.  Schematic  of  of  the  simultaneous  Velocity-Wavefront  measurements  using  Shack- 
Hartmann  2-D  wavefront  sensor  at  Notre  Dame. 
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The  experimental  set-up  at  Notre  Dame,  shown  in  Figure  5.12,  was  similar  to  the  one  at 
Caltech  with  two  major  differences.  Unlike  the  1-D  wavefronts,  measured  with  the  Malley  probe 
at  Caltech,  full  circular  2-D  wavefronts,  resolved  in  both  the  streamwise  (x)  and  the  spanwise  (z) 
directions,  were  collected  using  a  high-speed  Shack-Hartmann  wavefront  sensor.  The  wavefront 
aperture  was  50  mm.  The  second  difference  was  that  both  boundary  layers  on  the  opposite  tunnel 
walls  were  optically-aberrating.  In  order  to  minimize  the  optical  effect  of  the  upper  boundary 
layer,  a  Large  Eddy  Break  Up  (LEBU)  device  1  =  12  mm,  h  =  11  mm  device  was  mounted  on  the 
upper  wall  of  the  BL  development  section  to  reduce  the  level  of  optical  aberrations  caused  by  the 
upper  BL.  Its  trailing  edge  was  located  2.6  cm  upstream  of  the  start  of  the  optical  window. 
Estimations  based  on  results  from  Chapter  2  showed  that  the  LEBU  device  reduced  the  aero- 
optical  distortions  of  the  upper  boundary  layer  by  a  factor  of  1.7. 

Velocity  data  were  acquired  at  a  rate  of  2  kHz  using  a  commercially  available  2-D  PIV 
LaVision  system  in  double-pulse,  double-frame  mode.  The  laser  sheet  was  aligned  along  the 
tunnel  centerline,  and  the  frame  resolution  of  the  camera  was  768x768,  with  a  field  of  view  of 
approximately  100  mm.  The  time  interval  between  laser  pulses  for  each  measurement  was  10  ps. 
The  PIV  image  pairs  were  cropped  and  processed  in  DaVis  8.2  in  order  to  calculate  the  velocity 
vector  field  at  intervals  of  0.53  mm  in  the  streamwise  and  wall-normal  directions. 

Wavefront  measurements  with  the  spatial  resolution  of  40x40  subapertures  were 
simultaneously  acquired  using  a  Shack-Hartmann  sensor,  which  consisted  of  a  high-speed  camera 
(Phantom  vl611)  with  a  mounted  lenslet  array,  3  ps  after  the  first  laser  pulse,  so  as  to  have 
wavefront  data  that  closely  corresponded  to  each  image  pair  used  to  calculate  the  velocity  field. 
The  shutter  duration  for  the  Shack-Hartmann  sensor  was  only  0.452  ps,  and  the  points  at  which 
wavefront  measurements  were  obtained  were  distributed  at  1.2  mm  intervals  in  the  streamwise 
and  spanwise  directions.  The  aperture  diameter  of  the  wavefront  beam  was  approximately  5  cm, 
with  its  upstream  edge  passing  through  the  optical  window  at  x  =  133  cm. 

The  boundary  layer  thickness,  S99,  at  the  measurements  location  was  found  to  be 
approximately  15.6  mm.  Integration  of  the  velocity  profiles,  collected  with  PIV  system  gave  <5*  to 
be  2.4  mm  and  <9  =  1 .74  mm.  Based  on  these  values,  ReT  =  4,780  and  Re<?  =  15,500. 

Time-averaged  levels  of  OPD  from  the  velocity  data  were  found  approximately  25%  less 
than  the  OPDrms,  calculated  from  wavefronts.  The  difference  was  attributed  to  the  presence  of 
the  optically- weakened  upper  boundary  layer. 

5.8  Results:  Notre  Dame 

For  each  velocity  field,  the  estimated  1-D  slice  of  the  wavefront  was  calculated  using 
equations  (1.1)  and  (5.1).  The  resulted  wavefronts  (not  shown)  were  compared  to  the  1-D  slices  of 
the  measured  wavefronts  along  the  centerline  (Figures  5.13  and  5.14).  Instances  with  significant 
local  differences  between  the  estimated  and  actual  wavefronts  were  investigated.  As  the  estimated 
wavefronts  rely  on  SRA,  these  instances  should  correspond  to  the  cases  where  the  pressure/total 
temperature  fluctuations  are  not  negligible.  Two  events  in  the  measured  wavefront  slices  were 
considered,  local  OPD  minima  and  maxima,  and  the  velocity  fields  were  inspected  in  attempt  to 
identify  the  probable  relation  between  the  wavefront  minima/maxima  and  a  corresponding 
velocity  features. 

5.8.1  Local  OPD-minima. 

Several  representative  flow  fields,  when  the  local  wavefront  minimum  was  observed,  are 
shown  in  Figure  5.13.  The  upper  plot  shows  the  actual  2-D  OPD  (in  false  color)  and  the  1-D  OPD 
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slice  (black  line).  The  bottom  plot  shows  the  instantaneous  17-velocity  with  the  mean  convection 
velocity  of  optical  disturbance,  0.82t7x,  subtracted  (velocity  vectors),  superimposed  with  the 
instantaneous  fluctuating  u-field  (false  color).  It  was  observed  that  in  most  cases  the  local  OPD 
minimum  corresponds  to  the  presence  of  the  large-scale,  order  of  d,  structure,  which  is 
characterized  by  negative  w-fluctuations  near  the  wall  and  positive  ^-fluctuations  away  from  the 
wall.  Such  large-scale  structure  has  been  observed  to  correspond  to  the  signature  of  a  packet  of 
hairpin  vortical  structures  by  many  other  researchers  (Adrian,  2007,  Hutchins  et  al.,  2005,  for 
instance).  There  is  some  evidence  of  vortical  motion  reminiscent  of  hairpin  heads  in  the  PIV 
images. 
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Figure  5.13.  Top:  Representative  wavefronts  with  local  minima  near  the  PIV  centerline  (dashed 
line).  Measured  2-D  OPD  are  given  by  false  color  map,  blue  are  negative  values  and  red  are 
positive  values,  and  the  measured  1-D  OPD(x,0 )  slices  are  shown  as  black  lines.  Bottom:  the 
corresponding  velocity  fields.  Velocity  vectors  are  the  instantaneous  I/-velocity  with  0.82 /To 
subtracted,  superimposed  with  the  instantaneous  fluctuating  w-field  (false  color  map,  blue  are 
negative  values  and  red  are  positive  values).  Solid  magenta  line  indicates  a  v-location,  where  the 
local  mean  /7-velocity  is  0.82  of  the  freestream  speed. 


The  discrepancy  between  the  measured  and  actual  wavefronts  inside  these  large-scale 
structures  suggests  that  the  pressure  variations  are  not  small  inside  these  structures.  It  has  been 
shown  that  the  pressure  fluctuations  in  shear  layers  with  well-defined  vortical  structures 
significantly  contribute  to  the  overall  aero-optical  distortions  (Fitzgerald  and  Jumper,  2004).  So  it 
is  plausible  to  assume  that  a  packet  of  vortical  structures  in  the  boundary  layer  might  have  an 
associated  lower  pressure  region  inside,  leading  to  the  observed  discrepancy  between  the  SRA 
and  measured  OPD. 

As  wavefronts  are  integrals  of  the  density  field,  analysis  of  the  spanwise  extent  of  the 
local  minima  provides  non-intrusive  optical  measurements  of  the  instantaneous  spanwise  size  of 
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the  large-scale  structures.  For  instance,  from  Figure  5.13,  a  typical  spanwise  extent  is  found  to  be 
approximately  0.5fi?  with  a  typical  streamwise  extend  of  ~15.  Again,  these  observations  of 
predominantly  streamwise-elongated  large-scale  structures  are  consistent  with  findings  in  the 
literature  (e.g.  Hutchins  et  al,  2005). 


5.8.2  Local  OPD-Maxima. 


A  similar  analysis  was  performed  to  find  a  flow  structure  that  corresponds  to  the  local 
OPD-maximum.  Several  examples  of  velocity  fields  with  related  OPDs  are  shown  in  Figure  5.14. 
The  velocity  field  vectors  are  plotted  the  same  way,  as  in  Figure  5.13,  but  this  time  a  false-color 
map  shows  the  v-fluctuations.  It  was  observed  that  the  most  local  OPD  maxima  are  related  to 
compact  regions  of  large  positive  v-fluctuations,  residing  mostly  in  the  outer  portion  of  the 
boundary  layer.  One  plausible  explanation  of  why  these  ejection-like,  positive  vertical-velocity 
regions  corresponds  to  the  local  OPD  maximum  is  that  these  vertical  velocities  will  increase  the 
local  turbulent  thickness  of  the  boundary  layer  and  form  a  “bulge”  in  the  boundary  layer  interface; 
consequently,  the  upper  limit  of  integration  in  equation  (1.1)  and  therefore  the  OPD- value  would 
increase.  As  no  obvious  large-scale  vortical  structure  was  observed  in  these  representative 
velocity  fields,  local  OPD- maxima  are  most  probably  related  to  the  local  increase  of  the  boundary 
layer  thickness. 


Figure  5.14.  Top:  Representative  wavefronts  with  local  maxima  near  the  PIV  centerline  (dashed 
line).  Measured  2-D  OPD  are  given  by  false  color  map,  blue  are  negative  values  and  red  are 
positive  values,  and  the  measured  1-D  OPD(x, 0)  slices  are  shown  as  black  lines.  Bottom:  the 
corresponding  velocity  fields.  Velocity  vectors  are  the  instantaneous  17-velocity  with  0.8217* 
subtracted,  superimposed  with  the  instantaneous  fluctuating  v-field  (false  color  map,  blue  are 
negative  values  and  red  are  positive  values). 
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Again,  assuming  that  the  streamwise  and  the  spanwise  size  of  the  wavefronts  at  the  local 
maxima  are  related  to  the  size  of  these  vertical-velocity  regions,  wavefronts  in  Figure  6  suggest 
that  these  interface  “bulges”  are  approximately  0.7 S in  both  streamwise  and  spanwise  directions. 

Finally,  instances  where  the  measured  and  the  actual  wavefronts  agree  were  investigated 
(not  shown).  For  these  instances,  the  instantaneous  SRA  appears  to  be  valid  (pressure/total 
fluctuations  are  small)  and  no  significant  vortical  activity  in  the  boundary  layer  was  observed. 

5.9  Conclusions:  Notre  Dame 

A  complimentary  experimental  technique,  which  uses  simultaneous  velocity-wavefront 
measurements,  is  proposed.  Since  wavefronts  are  proportional  to  the  integrated  density  field,  it 
was  shown  that  the  comparison  between  the  velocity  field  and  the  wavefront  provide  additional 
information  about  the  large-scale  structures  inside  the  boundary  layer.  Plausible  relations  between 
different  structures  in  the  boundary  layer  and  the  corresponding  wavefronts  were  presented  and 
discussed.  The  SRA  was  used  to  investigate  the  relationship  between  optical  distortion  and  the 
velocity  field,  comparing  the  OPD  implied  by  the  velocity  field  with  the  one-dimensional  OPD 
from  the  Shack-Hartmann  sensor  at  the  location  of  the  PIV  plane.  2-D  wavefronts  were  shown  to 
provide  additional  information  about  the  streamwise/spanwise  statistics  of  the  large-scale  BL 
structures.  By  comparing  actual  wavefronts  with  the  estimated  ones,  using  the  Strong  Reynolds 
Analogy,  the  effect  of  the  local  pressure  field  and  the  underlying  vortical  structure  can  be 
estimated  and  studied. 
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